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Analysis  of  Electromagnetic  Interaction  with  Ships 
on  the  Ocean  Surface 

Jin  Au  Kong 

Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


1.  Research  Objective 

The  objective  of  this  research  project  has  been  to  develop  electromagnetic  models  for 
three-dimensional  (3-D)  objects  above  a  rough  surface,  analyze  and  predict  their  scattering 
responses,  determine  dominant  wave  interaction  mechanisms,  investigate  the  origins  of 
these  scattering  features,  and  provide  physical  interpretation  for  both  analysis  and 
synthesis  of  the  response. 


2.  Background 

Electromagnetic  scattering  from  surface  ships  is  an  important  subject  of  study  for  the 
Navy.  For  all  practical  purposes,  the  main  concerns  are  the  identification  from  radar 
signatures  and  reduction  of  radar  signatures  by  better  designs.  Both  of  these  require 
fundamental  understanding  of  how  electromagnetic  wave  interacts  with  the  ship  body  in 
the  presence  of  ocean  surface. 

In  the  past,  progress  in  this  area  has  been  hampered  by  the  complexity  of  the  problem.  A 
real  ship  is  composed  of  large  number  of  different  parts,  each  with  distinctive  response  in 
the  time  and  the  frequency  domains.  Moreover,  the  presence  of  rough  and  highly 
reflective  ocean  surface  creates  huge  amount  of  multiple  electromagnetic  interactions  with 
the  ship  body.  This  is  especially  true  at  low  elevation  angles,  where  most  practical 
applications  are  concerned.  Although  recent  advances  in  computational  electromagnetics 
and  increasing  capability  of  digital  computers  have  made  large  scale  electromagnetic 
scattering  problem  more  manageable  [1-3],  analytical  methods  [4-6]  are  still  needed  to 
facilitate  the  interpretation  of  numerical  solution.  Furthermore,  analytical  methods  can 
also  be  efficiently  combined  with  numerical  techniques  to  increase  the  overall 

computational  efficiency.  This  research  emphasizes  the  use  of  hybrid 

analytical/numerical  techniques  for  the  analysis  of  electromagnetic  interaction  with  surface 
ships. 


3.  Research  Achievements 

Numerical  simulation  techniques  for  electromagnetic  wave  scattering  by  arbitraiy  shaped 
objects  in  free  space  are  well  developed  using  wire  and  surface-patch  models.  The  theory 
and  numerical  approaches  associated  with  objects  near  flat  interfaces  of  layered  media 
have  also  been  studied  extensively  by  many  researchers.  However,  in  the  case  of  a  ship  on 


the  ocean  surface,  the  interface  is  inherently  rough  (sea  surface)  and  relatively  little  work 
has  been  reported.  In  principle,  the  standard  Method  of  Moments  (MoM)  can  be  used  to 
solve  for  the  unknowns  (surface  currents)  both  on  the  object  and  on  the  rough  surface. 
However,  the  need  to  discretize  the  rough  surface  significantly  increases  the 
computational  resource  requirements  compared  to  the  case  of  scattering  from  the  object 
alone. 

The  main  effort  in  our  work  has  been  the  development  of  iterative  techniques  that 
perform  perturbational  expansion  to  the  extended  boundary  condition  (EBC)  integral 
equations.  After  separating  the  equations  into  different  orders,  they  are  then  solved  using 
the  MoM  technique.  The  layered  dyadic  Green’s  function  formulation  is  applied  to  all 
orders  to  allow  only  unknown  current  distribution  on  the  discrete  conducting  object  but 
not  that  on  the  rough  surface  to  be  solved,  thus  saving  significant  amount  of  computer 
memory  and  CPU  time.  This  hybrid  SPM/MoM  overcomes  the  difficulties  of  having  to 
deal  with  a  large  randomly  rough  surface  that  the  standard  MoM  encounters  when  it  is 
applied  to  full  three-dimensional  problem.  The  zeroth  order  response  in  the  hybrid 
SPM/MoM  technique  has  been  found  to  be  equivalent  to  scattering  in  the  presence  of  a 
flat  interface,  and  the  first  and  higher  order  responses  can  be  solved  from  the  same  set  of 
(impedance)  matrix  equations  with  different  equivalent  sources  derived  at  the  ocean-air 
interface. 


3.1  Phase  I 

In  the  first  phase,  we  have  focused  on  the  development  of  far-field  interaction  models  for 
flat  plates  and  circular  cylindrical  shapes.  An  iterative  solution  for  scattering  from  a 
conducting  object  above  a  randomly  rough  dielectric  surface  has  been  formulated  using 
Huygens’  principle.  This  model  incorporates  analytical  techniques.  Small  Perturbation 
Method  (SPM)  and  Kirchhoff  Approximation  (KA)  for  rough  surface  scattering  and 
Physical  Optics  (PO)  for  scattering  from  the  conducting  body,  to  examine  the  direct  and 
multiple  scattering  mechanisms  for  the  coupled  configuration. 

We  first  considered  scattered  field  from  the  conducting  body  and  from  the  rough  surface 
alone.  The  secondary  scattering  components  were  evaluated  by  projecting  the  primary 
scattering  fields  to  the  discrete  object  and  rough  surface  as  additional  sources.  Using  the 
iterative  technique,  we  can  evaluate  the  spread  of  scattering  energy  associated  with 
different  surface  roughness.  The  main  lobe  in  bistatic  return  was  seen  to  decrease  with 
increasing  roughness,  whereas  some  sidelobes  pick  up  strengths.  This  spread  is  expected, 
as  the  scattered  energy  is  directed  away  from  the  specular  direction. 

A  numerical  technique  based  on  the  Method  of  Moments  (MoM)  solution  of  surface 
integral  equations  was  developed  to  validate  the  iterative  solution.  Monte  Carlo 
simulations  over  about  50  realizations  of  surface  profiles  were  used  to  evaluate  the 
statistical  average  return  in  the  combined  rough  surface/discrete  object  environment.  The 
surface  roughness  statistics  and  the  incident  waveform  were  assumed  to  be  the  same  as  in 
the  iterative  method.  Both  models  were  applied  to  investigate  the  effects  of  radar 
parameters  (incidence  and  observation  angles,  polarization  and  frequency),  and  target 


parameters  (size  and  location  of  the  object,  roughness  of  ocean  surface)  on  the  scattered 
fields.  It  was  observed  that  the  interactions  strongly  depend  on  the  frequency,  location  of 
object,  and  the  ocean  roughness.  Discrepancies  were  observed  between  the  iterative 
technique  using  Kirchhoff  approximation  and  MoM  when  the  object  is  slightly  above  the 
rough  surface  or  partially  submerged  (practical  case  of  interest),  which  suggest  that 
multiple  interactions,  beyond  the  secondary  scattering  covered  by  the  iterative  method, 
need  to  be  included.  The  use  of  the  hybrid  technique  outlined  in  the  next  sections 
improves  the  approximation  significantly.  As  a  by-product  of  this  study,  we  also 
examined  angular  correlation  function  of  scattered  waves  based  on  the  results  of  MoM.  It 
was  observed  that  the  detection  of  discrete  object  in  the  presence  of  rough  interface  would 
be  easier  if  angular  correlation  is  taken  along  certain  directions. 


3.2  Phase  II 

In  this  phase  of  research  [7,8],  we  have  developed  a  spatial-domain  formulation  for 
electromagnetic  scattering  from  a  rough  surface  for  an  arbitrarily  shaped  incident  wave 
using  the  extinction  theorem.  By  expanding  the  Green’s  function  and  surface  variables  in 
terms  of  the  rough  surface  height  function,  the  electric  field  integral  equations  (EFIE)  has 
been  decomposed  into  different  orders.  It  has  been  found  that  the  higher  order  equations 
are  of  the  same  form  as  the  zeroth  order  equations,  except  for  the  source  terms.  Therefore, 
for  each  order,  the  scattered  field  from  the  rough  surface  can  be  formulated  as  the  radiation 
fields  from  equivalent  sources  on  the  (flat)  mean  surface.  In  other  words,  we  only  need  to 
solve  the  zeroth  order  equations  with  different  sources  successively.  Moreover,  the 
problem  at  each  order  corresponds  to  a  flat  surface  problem,  which  can  be  solved  without 
discretizing  the  surface.  Up  to  the  first  order,  we  have  shown  that  the  total  field  of  the 
equivalent  source  on  a  PEC  rough  surface  is  the  field  scattered  by  the  rough  surface. 

The  advantage  of  using  the  spatial-domain  formulation  instead  of  the  frequency-domain 
formulation  as  in  the  conventional  SPM  technique  is  that  there  is  no  need  to  integrate  over 
the  spatial  frequency  .  Furthermore,  using  the  spatial-domain  formulation,  the  non- 
averaged  scattered  fields  are  well  defined.  However,  the  derivation  of  the  total  scattered 
field  of  higher  order  (n>2)  or  the  field  scattered  from  a  dielectric  rough  surface  requires 
more  analytical  effort  since  Sommerfeld  integrals  are  involved.  In  the  case  of  a  PEC  rough 
surface,  this  extra  analytical  effort  is  avoided  because  the  Green’s  function  can  be 
obtained  by  image  theory  and  does  not  involve  Sommerfeld  integrals.  By  considering  a 
plane  incident  wave,  we  have  shown  that  the  total  field  of  the  equivalent  source  is  the 
same  as  the  field  scattered  by  a  rough  surface  obtained  by  using  the  conventional  SPM 
technique. 

Particularly  for  a  rough  PEC  surface,  the  proposed  spatial-domain  technique  is  more 
powerful  than  the  conventional  SPM  for  the  calculation  of  the  scattered  field  because: 

(a)  The  incident  wave  is  not  restricted  to  be  a  plane  incident  wave  as  in  the  conventional 
SPM,  thereby  allowing  lower  grazing  angles  in  numerical  simulations. 

(b)  The  spatial  domain  formulation  can  be  used  to  calculate  the  near  field,  which  is 
particularly  important  to  account  for  near  field  interactions. 


(c)  The  non-averaged  complex  scattered  field  in  the  spatial-domain  formulation  does  not 
requires  the  evaluation  of  a  Fourier  integral  over  the  rough  surface,  therefore  is  a  well- 
defined  quantity. 

(d)  The  spatial-domain  formulation  is  more  computationally  efficient,  since  there  is  no 
integration  over  the  spatial  frequency  as  in  the  conventional  SPM. 

More  detailed  research  results  of  Phase  II  are  presented  in  the  Appendix  A  [8]. 


3.3  Phase  III 

In  this  phase  of  research  [7,9],  we  developed  a  hybrid  spatial-domain  SPM/MoM 
technique  to  calculate  the  EM  scattering  from  arbitrary  3-D  conducting  objects  above  a 
rough  surface  based  on  the  results  of  Phase  II. 

In  this  hybrid  technique,  the  Green’s  function  and  surface  variables  are  expanded  (using  a 
Taylor’s  expansion)  in  terms  of  the  surface  height  functions  on  the  (flat)  mean  surface. 
The  electric  field  integral  equations  based  on'  the  extinction  theorem  and  the  surface 
boundary  conditions  are  then  decomposed  into  the  different  orders.  Each  order  represents 
an  equivalent  flat-surface  scattering  problem  with  the  same  geometry  but  different 
equivalent  sources.  As  a  result,  the  problems  of  electromagnetic  scattering  from  an  object 
above  a  rough  surface  can  be  solved  efficiently  by  using  the  MoM  incorporating  the 
dyadic  Green’s  function  for  layered  media.  More  importantly,  only  the  surface  currents 
on  the  object  need  to  be  found  numerically  by  the  MoM. 

A  nice  picture  emerges  because  the  separation  of  the  solution  into  different  orders  helps 
the  identification  and  characterization  of  the  individual  interaction  terms  between  the 
object  and  the  rough  surface.  The  separation  of  the  return  field  into  the  sum  of  individual 
interaction  terms  allow  us  to  identify  the  coherent  and  incoherent  return  field,  and  thus 
better  characterize  the  rough  surface  effects  quantitatively. 

We  have  also  analyzed  the  case  of  an  object  partially  buried  in  a  PEC  rough  surface.  In 
comparison  to  the  case  of  an  object  above  a  rough  surface,  we  have  found  the  following 
differences  [7]: 

(a)  The  integration  area  of  the  equivalent  sources  on  the  mean  surface  is  changed  to 
accommodate  the  fact  that  part  of  the  (infinite)  flat  surface  is  covered  by  the 
interacting  object. 

(b)  The  are  extra  terms  in  the  expression  for  the  equivalent  sources  on  the  mean  surface. 
These  terms  are  related  to  the  correction  of  the  equivalent  source  due  to  the  change 
on  the  exposed  area  on  the  object  (part  of  it  being  buried)  and  to  the  modification  of 
the  equivalent  source  due  to  the  projected  area  of  the  rough  surface  on  the  mean 
surface. 

Several  numerical  results  have  been  obtained.  In  particular,  we  have  investigated  the 
interaction  of  horizontal  conducting  cylinders  above  a  rough  surface  with  Gaussian  power 
spectrum  and  validated  the  results  against  standard  MoM  results.  Monte  Carlo  results 


have  also  been  generated  using  1 00  rough  surface  realizations.  The  number  of  unknowns 
for  typical  cases  considered  has  been  reduced  by  roughly  20  times.  Moreover,  a  full 
Monte  Carlo  monostatic  RCS  calculation  which  would  require  3  months  in  a  our  500 
MHz  Digital  Alpha  21164A  processor  using  the  standard  MoM,  was  reduced  to 
approximately  1 7  hours  using  the  hybrid  MoM/SMP  technique. 

More  detailed  research  results  of  Phase  III  are  presented  in  the  Appendix  B  [9]. 


3.4  Phase  IV 

A  basic  limitation  of  usual  numerical  techniques  to  analyze  the  problem  of 
electromagnetic  scattering  from  rough  surfaces  is  the  need  to  truncate  the  rough  surface 
because  of  limited  computational  resources.  This  leads  to  erroneous  results  due  to 
artificial  edge  truncation  when  ideal  plane  waves  are  used  to  excite  the  configuration.  The 
tapered  wave  concept  is  based  on  providing  an  illumination  for  the  numerical  simulation 
which  resembles  the  plane  wave  case  to  be  modeled  closely  to  the  center  of  the  scattering 
scenario  (including  the  choice  of  a  particular  p'olarization)  while  its  intensity  becomes 
negligibly  small  upon  approaching  the  artificially  introduced  edges  of  the  rough  surface. 
Therefore,  unwanted  edge  effects  due  to  the  primary  incident  wave  are  minimized  and  a 
proper  normalization  of  the  computed  scattering  coefficients  allows  a  meaningful 
comparison  with  the  ideal  plane  wave  case.  Moreover,  near  field  quantities  such  as 
current  distributions  induced  near  the  center  of  the  tapered  wave  are  also  expected  to  be 
similar. 

The  tapered  wave  should  be  constructed  in  such  a  way  that  it  satisfies  the  Maxwell’s 
equations.  This  helps  to  increase  the  confidence  in  the  results  obtained  from  numerical 
simulations.  Furthermore,  it  should  also  be  possible  to  substitute  it  for  a  plane  wave  of 
arbitrary  polar  and  azymuthal  angles  of  incidence  without  loss  of  polarization  or 
degradation  of  tapering. 

In  this  research  phase  [10],  we  have  considered  the  problem  of  constructing  a  3-D  tapered 
wave  with  dominant  polarization  state  for  all  angles  of  incidence  as  a  superposition  of 
plane  waves,  taking  into  account  both  propagating  and  evanescent  waves.  The  use  of  the 
simple  Gaussian  plane  wave  spectrum  was  considered  in  order  to  avoid  problems  near 
grazing  incidence.  Moreover,  we  have  proposed  a  choice  of  unit  vectors  which  leads  to 
integrand  which  are  continuous  at  normal  incidence.  As  a  result,  a  clean  footprint  (no 
degradation  of  tapering)  with  clear  polarization  (no  loss  of  dominant  polarization  state) 
has  been  obtained  for  all  angles  of  incidence.  Besides  removing  the  singularity  at  the  origin 
of  the  wavenumber  space,  the  proposed  special  choice  of  polarization  vectors  has  lead  to 
a  least  square  error  property  of  the  wave.  The  choice  of  polarization  vectors  was  shown 
to  lead  to  an  exact  solution  of  the  Maxwell’s  equations  which  is  an  optimal  approximation 
of  an  ideal  but  non-Maxwellian  tapered  field  that  is  constructed  by  multiplying  a  scalar 
tapered  wave  with  a  constant  polarization  vector.  Approximate  3-D  scalar  and  vector 
tapered  waves,  which  can  be  evaluated  without  any  numerical  integrations,  has  also  been 
derived. 


More  detailed  research  results  of  Phase  IV  are  presented  in  the  Appendix  C  [10]. 


4.  General  Conclusions  and  Future  Directions 

4.1  General  Remarks 

In  our  approach  we  successfully  combined  an  analytical  preprocessing  (perturbational 
expansion)  of  the  problem  of  electromagnetic  scattering  from  an  object  in  presence  of 
rough  surface  to  be  used  in  a  numerical  solution  by  the  MoM.  This  combination  of 
analytical  and  numerical  treatment  lead  not  only  an  improvement  on  the  overall 
computational  efficiency  of  the  resultant  (hybrid)  scheme,  but  also  permits  a  better 
physical  interpretation  of  the  results  given  by  the  computer  (i.e.,  separation  of  the  object- 
rough  surface  interaction  into  different  orders)  which  is  much  difficult  using  a  strictly 
numerical  procedure. 

The  development  of  theoretical  framework  needs  to  be  coupled  with  experiments,  both 
for  validation,  and  for  evaluating  new  measurement  techniques.  In  essence,  one  can 
perform  theoretical  simulations  that  take  into  account  the  measurement  setup  and  assess 
the  impact  of  testing  probes.  One  of  the  long-term  goals  of  our  research  is  to  provide 
recommendation  on  experimental  schemes. 

However,  many  important  issues  are  unresolved  still.  For  instance,  being  a  perturbational 
expansion,  the  hybrid  scheme  can  become  very  involved  for  non-slightly  rough  surfaces. 
Hybridization  with  other  techniques  and  use  of  composite  surface  models  need  to  be 
investigated  in  connection  with  more  general  rough  surface  conditions.  Moreover,  in  the 
course  of  our  research  we  have  dealt  only  with  the  forward  problem.  For  the  identification 
of  radar  signatures  and  reduction  of  radar  signatures  the  forward  techniques  need  to  be 
coupled  to  inverse  scattering  techniques.  These  and  other  issues  are  discussed  below. 

4.2  Future  Work 

One  immediate  future  goal  is  to  extend  the  applicable  range  of  our  hybrid  technique. 
Major  technical  issues  to  be  addressed  include  the  extension  of  the  hybrid  technique  to 
composite  rough  surface  model  with  impedance  boundary  condition,  which  will  be  used 
to  characterize  different  sea  states,  and  the  integration  with  standard  high-frequency 
techniques,  such  as  ray-optical  methods,  for  large  objects. 

4.2.1  Impedance  boundary  condition. 

To  obtain  a  closed-form  Green’s  function,  thus  avoiding  the  need  of  extra  integrations,  we 
have  assumed  that  the  ocean  surface  behaves  like  a  perfectly  conductor.  In  reality,  the 
behavior  of  highly  conducting  interface  is  quite  different  for  low-grazing  incidence  angles. 
Major  difference  is  exhibited  in  the  angular  variation  of  reflection  coefficient  for  the 
vertical  polarized  waves.  The  reflection  coefficient  converges  to  minus  one  (-1)  at  near 
grazing  incidence  even  at  the  limit  of  conductivity  reaching  the  infinity.  On  the  other 
hand,  the  reflection  coefficient  at  a  perfectly  conducting  interface  is  always  plus  one  (+1). 
This  non-uniform  convergence  property  can  conceivably  contribute  to  differences 


between  the  perturbational  series  expansion  for  the  perfect-conducting  interface  case  and 
that  for  the  high-conductivity  interface  (real  sea  water)  case.  Because  of  this  concern,  we 
propose  to  extend  our  formulation  to  incorporate  impedance  boundary  conditions  for 
highly  (non-perfect)  conducting  interfaces.  Doing  so  should  represent  the  sea  surface 
property  better  and  avoid  the  non-uniform  convergence  problem. 

4.2.2  Composite  rough  surface  model. 

Based  on  practical  experience,  the  effect  of  rough  ocean  surface  on  scattering  from  ship 
objects  becomes  non-negligible  only  at  high  frequencies.  In  such  cases,  it  is  necessary  to 
consider  a  two-scale  composite  rough  surface  representation  of  the  ocean  surface  [11]. 
The  perturbational  formulation  developed  here  can  be  easily  generalized  to  deal  with  the 
two-scale  model.  With  the  two-scale  model,  the  large  scale  underlying  surface  structure  is 
sufficiently  smooth  to  apply  high-frequency  scattering  techniques  such  as  physical  optics 
or  ray-optical  techniques.  The  ship  objects  can  also  be  facetized  and  included  in  the 
computational  domain  with  the  large-scale  rough  surface.  The  hybrid  technique  developed 
here  may  then  be  used  in  parts  of  the  object  and  surface  where  there  strong  multiple 
interactions  are  expected  to  improve  the  accuracy,  with  the  perturbational  formulation 
being  introduced  to  account  for  the  effect  of  small-scale  roughness. 

4.2.3  Finite-distance  source  synthesis. 

At  low  grazing  incidence,  it  is  important  to  observe  that  the  scattered  field  from  the  rough 
ocean  surface  will  interfere  and  merge  with  the  incident  field  that  is  propagating  near  the 
surface.  The  net  result  is  that  the  actual  incident  field  reaching  ship  objects  is  modified; 
likewise  the  return  field  to  the  sensor  will  be  modified  owing  to  scattering  from  the  rough 
surface.  To  deal  with  this  practical  problem,  vye  propose  to  include  a  source  term  at  a 
finite  distance  into  the  iterative  formulation  instead  of  a  synthesized  wave  beam  that  may 
not  represent  the  true  incident  field  near  the  ship  objects.  The  source  term  should  be 
constructed  to  emulate  the  transmitting  antenna  properties.  Taking  into  consideration 
the  statistical  properties  of  rough  ocean  surface  and  the  antenna  beam  pattern,  we  can 
calculate  the  surface  current  distribution  for  near-field  illumination  as  well  as  the 
dispersion  characteristics.  The  calculated  results  can  then  be  compared  against 
experiments  that  involve  objects  on  a  real  ocean  surface. 

4.2.4  Inverse  problem. 

The  essence  of  identifing  radar  signatures  can  be  considered  an  inverse  scattering  problem. 
The  work  here  begin  with  the  forward  scattering  model  developed  and  proceed  to  the 
inverse  model.  The  link  between  the  two  is  established  through  optimization  techniques. 
In  the  optimization  approach,  the  developed  forward  model  and  observed  measurements 
are  used  to  retrieve  unknown  source  parameters.  This  is  usually  the  most  appropriate 
inversion  method  when  the  direct  scattering  theory  is  very  complicated.  In  addition,  this 
method  offers  a  great  flexibility  in  the  choice  of  direct  scattering  models,  the  parameters  to 
be  inverted,  and  data  to  be  employed.  We  have  previously  applied  the  optimization 
technique  to  the  inversion  of  ground  truth  information  from  satellite  remote  sensing  data, 
such  as  inversion  of  sea  ice  parameters  from  the  CRRELX  experiment  [12].  One 
drawback  of  the  optimization  technique  is  the  amount  of  computational  resources 
required,  since  multiple  forward  scattering  solutions  are  needed.  One  way  to  overcome 
this  problem  is  to  extract  reduced-order  models  from  the  complete  forward  scattering 


solution  and  use  them  at  intermediate  stages  of  the  optimization  procedure.  An  example 
of  reduce-order  modeling  consists  in  the  replacement  of  the  scattering  object  by  scattering 
centers  [13]. 

Aside  from  the  optimization  techniques,  direct  time-domain  inversion  techniques  [14,15] 
can  also  be  considered.  For  example,  in  the  singularity  expansion  method  (SEM)  [16] 
feature  extraction  is  achieved  by  studying  the  natural  frequencies  contained  in  the  time- 
domain  response.  This  technique  might  well  be  suitable  for  identifying  local  dispersion 
characteristics  of  substructure  components. 

For  the  inverse  scattering  problem,  non-uniqueness  and  inversion  stability  issues  usually 
arises  with  noisy  data.  It  is  possible  to  reduce  the  effects  of  these  uncertainties  by  using 
multi-dimensional  data  like  multi-frequency,  multi-angle,  and  polarimetric  data.  For 
instance,  electromagnetic  wave  polarization  returns  are  particularly  sensitive  to  the 
geometry  and  the  degree  of  heterogeneity  of  scattering  objects.  Using  polarimetry,  the 
copolarized  and  depolarized  returns  for  an  arbitrary  polarized  incident  wave  can  be 
measured  and  exploited  for  the  classification  and  identification  of  targets  [17]. 
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Abstract 

We  present  a  derivation  of  a  spatial-domain  formulation  for  electromagnetic  wave  scattering  from 
a  rough  surface  for  an  arbitrary-shaped  incident  wave  using  the  extinction  theorem.  By  expanding  the 
Green’s  function  and  surface  variables  in  terms  of  the  rough  surface  height  function,  the  electric  field 
integral  equations  (EFIE)  are  decomposed  into  different  orders.  For  each  order,  the  scattered  field  of 
rough  surface  is  formulated  as  the  radiation  of  equivalent  sources  on  the  mean  surface.  The  spatial-domain 
formulation  is  simpler  and  not  restricted  to  use  a  plane  incident  wave  as  conventional  small  perturbation 
method  (SPM)  does.  The  spatial-domain  formulation  can  also  provides  a  non-averaged  scattered  field 
without  using  the  Fourier  transform  of  the  rough  surface.  The  numerical  result  for  a  non-planar  incident 
wave  is  verified  using  the  method  of  moments  (MoM). 

Keywords 

Rough  surface,  electromagnetic  wave  scattering,  spatial  domain,  equivalent  source,  small  perturbation 
method. 


I.  Introduction 

The  electromagnetic  wave  scattering  from  a  slightly  rough  surface  has  been  studied 
extensively  using  (1)  the  Rayleigh  hypothesis  [1],  [2],  and  (2)  the  extinction  theorem  [3], 
[4].  In  all  of  these  methods,  a  plane  incident  wave  has  been  assumed  so  that  the  scattered 
fields  can  be  obtained  by  expansion  in  terms  of  plane  waves.  Therefore  the  scattered 
fields  are  the  superposition  of  plane  waves.  The  averaged  scattered  fields  are  calculated 
by  ensemble  average  over  the  rough  surface,  thus  the  expression  after  the  average  is  the 
integral  over  the  wavenumber  for  a  function  of  the  power  spectrum  density  of  the  rough 
surface.  The  limitation  of  the  conventional  SPM  is  obvious:  (1)  a  plane  incident  wave  must 
be  assumed,  which  may  arise  the  difficulty  when  a  incident  wave  with  a  grazing  angle  is 
considered,  (2)  a  Fourier  integral  of  a  rough  surface  appears  in  the  non-averaged  scattered 
field,  and  (3)  it  is  time  consuming  to  evaluate  the  integration  over  the  wavenumber  when 
the  closed  form  of  the  scattered  field  could  not  be  obtained. 

In  this  paper,  we  derive  the  formulation  in  spatial  domain  for  the  field  scattered  by  a 
slightly  rough  surface.  In  Section  II,  based  on  the  extinction  theorem,  we  formulate  the 
equivalent  sources  on  the  mean  surface  for  a  dielectric  rough  surface.  Then  we  introduce 
the  iterative  scheme  to  calculate  the  total  scattered  field.  In  Section  III,  in  order  to 
calculate  the  reflected  and  transmitted  fields  of  the  equivalent  sources  on  the  mean  surface, 
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we  expanded  the  Green’s  function  as  the  integral  in  frequency  domain.  The  closed  form 
in  spatial  domain  of  the  first-order  field  scattered  by  a  perfectly  conducting  (PEC)  rough 
surface  is  derived  in  the  end  of  this  section.  In  Section  IV,  we  show  that  the  first-order 
scattered  field  in  spatial  domain  from  a  PEC  rough  surface  is  the  same  as  the  conventional 
SPM  by  assuming  a  plane  incident  wave.  For  arbitrary-shaped  incident  wave  such  as 
Gaussian-tapered  wave,  the  spatial-domain  formulation  provides  accurate  result  which  is 
verified  using  the  method  of  moments  (MoM)  in  Section  V. 

II.  Electromagnetic  Field  Integral  Equations 

In  the  recent  investigation  on  the  electromagnetic  wave  interaction  between  a  perfectly 
electric  conducting  (PEC)  object  and  a  rough  surface,  it  was  found  that  the  problem  can 
be  solved  by  replacing  the  rough  surface  with  a  flat  one  and  placing  an  equivalent  surface 
on  the  fiat  mean  surface  [5],  [6].  It  immediately  arise  the  question  that  if  the  equivalent 
source  method  can  be  applied  to  calculate  the  electromagnetic  wave  scattered  by  a  plain 
rough  surface. 

Consider  an  electromagnetic  wave  E{(r)  incident  upon  a  rough  surface  Sr  with  surface 
height  function  z  =  f(f'±)  and  mean  surface  S0  coincident  with  the  xy-plane  as  illustrated 
in  Fig.  1.  The  upper  and  lower  spaces  \\  and  Vi  are  homogeneous,  isotropic  media 
characterized  by  (ei, /q)  and  (e2,//2),  respectively.  By  applying  the  extinction  theorem, 
the  electric  field  integral  equations  are  written  as,  respectively,  for  the  observation  point 
at  f  €  V2 

E*(f)  +  f  dS'  {iuviiGi(f,  f7)  •  [n^r7)  x  /^(f7)] 

J  Sr 

+V  x  Sftf')  •  [^(f7)  x  ^(f7)]}  =  0  (1) 

and  for  r  G  V2 

Js  dS'  {iu>H2G2(f,f')  •  [^(f7)  x  /f2(f7)] 

+V  x  G2(f, f7)  •  [n2(fO  x  E2(f7)] }  =  0  (2) 

where  G\  and  G2  are  dyadic  Green’s  functions  for  unbounded  regions  Vi  and  V2,  respec¬ 
tively.  ni(f')  denotes  the  local  normal  vector  pointing  from  the  rough  surface  Sr  to  the 
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Fig.  1.  Electromagnetic  wave  scattering  by  a  rough  surface. 

upper  region^,.  The  normal  vector  n^)  is  anti-parallel  to  *,(?),  i.e., 

f  l(7  “d  flM  e'eCtric  “d  ma«"dic  fid*  on  the  rough  surface  in  region  V„  while 
E2(r )  and  H-Jr )  axe  surface  fields  on  Sr  in  region  V2.  If  the  region  V2  is  dielectric,  then 
the  tangential  fields  are  continuous^  thus 

ni(f')  x  ff,( Y)  =  (f')  x  U2(?)  =  l^a(r'x)  (3) 

"i(f')  x  E^f'j  =  Ai(f')  x  E2(f')  =  ~b(fx)  (4) 

where  a(?x)  and  b(fx)  are  two  new  surface  variables  defined  on  the  mean  surface  S,  is 
the  intrinsic  impedance  of  the  upper  region  V„  i.e.,  *  =  ^  „  the  projMtio„ 

of  the  infinitesimal  area  dS'  on  the  mean  surface  With  the  new  surface  variables,  the 
integral  equations  are  rewritten  as 

Ei(r)  + df1 l  {f,  r')  •  a(f^) 

+V  x  f')  •  6(f^)}  =  0  for  f  €  H  (5) 

fSg  ^r-L  ^2~G2(f,  f')  •  a(f'M) 

+V  x  G2(f,  f')  •  6(fl)}  =0  for  f  e  Vi  .  (6) 

Theoretically,  given  a  deterministic  rough  surface  profile  Sr,  the  unknown  surface  variables 
o(rx)  and  b(7x)  can  be  solved  from  Eqs.  (5)  and  (6).  For  the  special  case  in  which  the 
surface  is  fiat,  the  surface  variables  d(fx)  and  b(7x)  have  only  the  horizontal  components 
and  the  local  coordinate  f  in  the  dyadic  Green's  functions  can  be  replaced  by  f 
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A.  Expansion  of  Green’s  Function  and  Surface  Variables 

As  Illustrated  in  Fig.  2,  the  local  position  vector  ?  on  the  rough  surface  can  be  expressed 
as  the  sum  of  the  horizontal  vector  and  the  vertical  vector  zf(r'±).  Therefore  the  scalar 
Green  s  function  ga  in  region  Va  (where  a  =  1  or  2)  can  be  expanded  in  terms  of  the 


surface  height  function  f(r'±)  on  the  mean  surface  S0, 


tf(K) 


Fig.  2.  The  local  position  vector  r'  on  a  rough  surface  is  the  sum  of  the  horizontal  and  vertical  components. 


9a(r,r')  =  ga(r,r'±  +  zf(r'±)) 

OO  1  am 

=  (7) 

where  the  following  property  of  the  scalar  Green’s  function  has  been  used: 

dm  ft™ 

0^9a(f,  f)  =  (-l)m  g-^9a(f,  f1).  (8) 

Thus  the  dyadic  Green’s  function  can  be  expressed  as 

^o(f.r')  =  (7  +  -*fw)  gjf,  f'j 

=  (9) 

The  surface  variables  are  also  written  as  the  series  expansion  as  follows: 

OO 

aft)  =  E  [af’ftJ  +  ia^ft)]  (10) 

m= 0 
oo 

*ft)  =  E  [6im)ft)+4m,ft)l-  (ii) 

m=0 

In  Eqs.  (10)  and  (11),  the  separation  of  the  z-components  allows  the  only  unknown  surface 
variables  to  be  the  tangential  component  since  the  z-components  of  the  m-th  order  can 
be  expressed  in  terms  of  the  tangential  components  of  order  (m  —  1),  as  we  will  see  below. 
For  the  definition  of  a  and  b  as  in  Eqs.  (3)  and  (4),  the  following  identities  hold: 


(12) 

(13) 
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where  the  normal  vector  is  related  to  the  surface  height  function  by 

nil 


"l(rJ  “  |-Vi/(fl)+I|- 


(14) 


Substituting  the  series  expansion  of  a  and  b  as  in  Eqs.  (10)  and  (11)  into  Eqs.  (12)  and 
(13),  respectively,  and  assuming  that  the  derivatives  df(f'±)/dx'  and  df(f'±)/dy'  are  of  the 
same  order  as  fci/(fj_),  we  then  get  the  m-th  order  ^-components  of  the  surface  variables 
a  and  b  in  terms  of  their  (m  —  l)-th  order  tangential  components  as  follows: 

sf’ft)  =  Vl/(f'x).ai”-,>(fi)  (15) 

(16) 

B.  The  n-th  Order  Equation  and  Equivalent  Sources 


Substituting  the  series  expansion  of  the  dyadic  Green’s  function  [Eq.  (9)]  and  the  surface 
variables  [Eqs.  (10)  and  (11)]  into  the  integral  Eqs.  (5)  and  (6),  and  re-grouping  the  terms 
with  the  same  order  of  the  surface  height  function  we  obtain  the  following  n-th 

order  equations: 

MnHf)  +  f  dfj.  f'±)  ‘  Oj^l) 

J  So 

+ V  x  V, (f,  f'±)  ■  b™ (f'x)}=  0  for  f  6  V2  (17) 

+  /  dfx  { ik2—G2(f,  f^)  •  af  )(f(L) 

J  S0  ^  7)i 

+ V  x  V2  (F,  f'x)  •  >  (rl) }  =  0  for  f  e  Vi  (18) 

where  for  the  zero-th  order  (n  =  0),  5j0)(f)  =  Ei(f)  and  E^ (f)  =  0.  For  higher  orders 
(n  >  1),  the  first  terms  in  the  equations  above  are  written  as 

m=l  1 
pxm  _ 

gpj  {><:i3i(f,f'j.)  •  aj"m)(f'x) 

+  Vx51(f,f'x)-Bi“-™>(fi)} 

{<*,  G,(f,f'x)  ■  z  [V'x/( r-x)  ■  ai"-’"l(f'x)] 

+  v  x  5,  (f,  r-j  •  *  [v'j(?o  ■  bfrm\f'0] }  (19) 
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and 

®S?(*0  =t±f  <  [-/(r-l)]"1 

m=l  m!  7 

+  vx52(f,fi)4"-"Vi)} 

^  {<*4 1  G2(f ,  fl)  ■  z  [V'x/(fi)  ■ 

+  V  x  ff2(f ,  fi)  •  I  [V'x/(fl)  •  Jf->(fU] }  .  (20) 

Comparing  the  zeroth  order  (n  =  0)  integral  equations  with  Eqs.  (5)  and  (6)  for  the 
flat  interface  [/(fjJ  =  0],  we  notice  that  they  are  the  same  except  for  an  additional 
superscript  (0)  at  the  surface  variables  a±  and  b±.  Therefore  the  solutions  and  for 
the  zeroth  order  equations  should  be  the  same  as  Eqs.  (5)  and  (6)  for  which  the  surface 
is  considered  to  be  flat.  For  the  higher  order  equations  (n  >  1),  we  also  notice  that  they 
are  in  the  same  form  as  the  equations  for  flat  surface,  except  for  the  substitution  of  £* 
by  E\n\  and  the  additional  “source”  term  in  the  lower  region  V2.  We  call  the  terms 
E\n^  and  e\^  are  fhe  equivalent  sources  (ES)  placing  on  the  mean  surface  Sa  of  the  rough 
surface  as  illustrated  in  Fig.  3.  By  denoting  the  corresponding  reflected  field  as  for 
the  equivalent  source  E in  the  upper  region  V\ ,  and  the  transmitted  field  as  E^  for  the 
-equivalent  source  ~E^  in  the  lower  region,  the  total  scattered  field  from  the  rough  surface 
illuminated  by  the  incident  wave  E{  can  be  expressed  by 

e,(t)  =  £  (ifV) + eJ’V))  +  E  ES’(f).  (2i) 

n=l  V  '  n=l 


C.  Iterative  Scheme  for  Total  Scattered  Field 

Eqs.  (19)  and  (20)  provide  the  spatial-domain  formulation  of  the  electric  field  from 
the  equivalent  source  in  terms  of  the  lower  order  surface  variables.  The  following  is  the 
iterative  procedure  to  calculate  the  scattering  field  from  the  rough  surface: 
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Fig.  3.  Radiation,  reflection  and  transmission  of  equivalent  sources  on  the  mean  surface. 


1.  The  zero-th  order  scattered  field:  Calculate  the  reflected  ET(f)  and  Hr(f)  of  the  incident 

fields  Ei(r)  and  #,-(r)  from  the  mean  surface.  Thus  the  zero-th  order  scattered  field  is  the 
same  as  the  reflected  field  E^\f)  =  Er(f)  and  H^\f)  =  Hr(f). 

2.  The  first-order  scattered  field:  Calculate  the  total  zero-th  order  fields  on  the  mean 
surface 


E<0)(f)  =  Et(  r)+B<,0)(f) 

b <0)w  =  H((f)+^o)(f). 


(22) 

(23) 


(24) 

(25) 


Then  calculate  the  surface  variables  of  the  zero-th  order  can  be  obtained  by 

4Vi)  =  t x  if'Vi) 

4Vx)  =  ixS"Vj. 

Thus  the  first-order  scattered  field  can  be  calculated  by 

+  (26) 

sf’tf)  =  B<,)(f)+B<r,,(f)  +  5<«(f)  (27) 

where  Bi  (f)  and  if,  \f)  are  the  fields  of  the  first-order  equivalent  source  which  can  be 
calculated  by  using  Eq.  (19)  for  n  =  1.  E^f)  and  #>)  are  the  reflected  fields  of  the 

equivalent  source  in  the  upper  region  Vi,  and  Bg’(f)  and  lf"(r)  are  the  transmitted  fields 

from  the  equivalent  source  in  the  lower  region  V2. 

3.  Calculate  the  n-th  order  scattered  field:  Calculate  the  surface  variables 

s±,(fx)  =Viz 


(28) 
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Thus  the  n-th  order  scattered  field  can  be  calculated  by 


=  Eli'\f)+E(;)  (f)+E!;’(f)  (so) 

+  +  (31) 


where  and  are  the  fields  of  the  n-th  order  equivalent  source,  and 

H^\f)  are  the  reflected  fields  of  the  equivalent  source  in  the  upper  region  Vi,  and 
and  H^\r)  are  the  transmitted  fields  from  the  equivalent  source  in  the  lower  region  V2. 


III.  Integral  Representation  of  Equivalent  Source 

The  integral  representation  of  dyadic  Green’s  function  in  region  Vj.  and  V2  can  be  written 
as  [7],  respectively, 

7^  ,  i 

Gi(r,r)  =  -zz  ,  0  + 


and 


k\  '  87T2 

J  dkj_^-Ai(ku)etkl<r~r,) 

[  dkxT-Sfii(-kiz)e**rif~r 

J  K\z 


z>  z' 
z  <  z' 


^  9Ar’f ')  i 

u2(r,  r  )  =  -22 — To —  + 


k\  '  87T2 

J  Koz 


?) 


z>  z' 


r 

1 1  d~k^^(-ki‘)e>K’<r 


f~7)  z  <  z< 


(32) 


(33) 


where  we  define 


•^i(^iz)  —  £i(fciz)ej(A;lz)  +  hi(kiz)hi(kiz)  (34) 

A2(k2z)  =  &‘i{k2z)e2{k2z)  +  h2{k2z)h2(k2z)  (35) 

and  ei)2  and  hl>2  are  unit  polarization  vectors  defined  in  Appendix.  The  up-going  wave 
vectors  kY  and  k2  are  expressed  as  k\  =  xkx  +  yky  +  zklz  and  k2  =  xkx  -1-  ykv  +  zk2z, 
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respectively,  and  the  down-going  wave  vectors  K\  and  K2  are  expressed  as  K\  =  xkx  4- 
ykv  —  zku  and  K2  =  xkx  +  yky  —  zk2z,  respectively.  Defining  the  dyads 

Bi{kiz)  =  kixAi(ku) 

=  -hi(kiz)ei(kiz)  +  ei(kiz)hi(klz)  (36) 

B2{k2z)  =  M  A2{k2z) 

=  —  h2(k2z)^2(k2z)  +  ^2{k2z)h2(k2z)  (37) 


we  write  the  expressions  for  V  x  Git2(r,r')  in  terms  of  the  integral  of  B\^.  The  operator 
dm/dzm  applying  on  Gii2(f,  f)  and  V  x  Git2(f,f')  creates  the  following  eigenvalues: 


gm 

dzm 


(iki,2z)m  .  z>  z' 

(~iki,2z)m  z  <  z'. 


(38) 


Therefore,  the  radiation  field  from  the  equivalent  source  on  top  of  the  mean  surface  can 
be  written  as  follows: 

For  up-going  wave  (z  >  z') 

m=l 

{/ dk±(iklzr~‘  ■  a'rm,( fl) 

+ / dkL  B,(^y -  s'rm,(fi)} 

[jdhiik.r-2  %(klz)^<'-<"> 

+ j dk^ih.f1-2 

■ilWJ't'rte])  (39) 


{-  [ dk±(ik2z)m-1  l2(k2z)eik><^  -atm)(r±) 
l*7i  J 
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+  J dkx(ik2,)m~' 

df'A~m)r' 
J ikxiih. 2l)m-2  52(Me‘M,-'') 

+  j  dk±  % (k2x)e'l‘  l,-,l} 


(40) 


For  down-going  wave  (z  <  z') 

!fdkJ.(-ik„)m-,em-('-r> 

+ f  dk±  (~iklz)m-‘  e®  J 

-  'gt^/rA-nur1 

•i[v'x/(fi)-ar”‘»(f'x)J 
+  J  dkx  (-ik,x)"‘-2  I,  (-ku)^'v-r] 

.|[v'x/(f'x).6i"'m>(f'x)]J  (41) 


m=l 

+  J  dk±  (~ik2z)m~l 
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B2(-k7,)  b'"  m| (f±) J 

f  dk±(-,k2zy-*  %(-k2,)eni^-n 

+  /  dkj.{-ih,r~2  ll(-fc«)e‘F*'(''"'', 

^[V'x/(f'J.)-5f-m»(f'x)]J.  (42) 

The  down-going  wave  and  up-going  wave  e\^  wiH  be  reflected  by  and  transmitted 
through  the  mean  surface,  respectively. 

A.  Reflected  Field  of  Equivalent  Source 

To  calculate  the  reflected  wave  of  the  equivalent  source,  we  decompose  E ^  into  a  sum 
of  TE  and  TM  waves.  Recalling  the  definition  for  Ai,  B i,  A2,  and  B2,  we  know  that  the 
product  of  ei(—kiz)ei(—kiz)  or  ei(— kiz)hi(— fclz)  with  any  polarization  will  be  a  TE  wave, 
and  the  product  of  hi(—klz)hi(—kiz)  or  hi(—klz)ei(—kiz)  with  any  polarization  will  be  a 
TM  wave.  To  get  the  reflected  wave  for  the  down-going  wave  for  z  <  zf,  we  do  the 
following: 

1.  Multiply  the  TE  wave  terms  by  RTE(kiz)  and  the  TM  wave  terms  by  R™(kiz)  in  the 
integral  representation  of  the  down-going  wave 

2.  Change  the  sign  of  Aqz  for  the  first  unit  polarization  vector  in  the  tensors  A\  and  Bu 
and 

3.  Change  the  down-going  wave  vector  Kx  to  the  up-going  wave  vector  ki  in  the  exponen¬ 
tial  terms. 

Thus  we  find 

{/ dk±  (-4*1.)”-’  (fire(*i,)«l(*i.)«i(-*u) 
+R™(klz)hl(klz)h1(-klz))  •  a(xn"m)(fl) 
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+  J  dk±  (~iku)m-'  (-fl™(»!u)Si  (*„)«,(-*„) 
+RTE(klz)el(klz)hl(-ku))  ■  6in-m)(fl)  J 

U  dkx(-rk„)’n-2 

(i?TE(fciz)ei  (fciz)ei  (-ku) 

+/2rw(fcxz)  Ai  (&iz)  hi  (— &iz)) 

•z[vl/(fl)-ain-m)(fl)] 

+  J  dk±  (—ik\z)m~2 

+i2rs(A:lz)e1(A:lz)Ai(— fciz)j 


■s  [v'x/(fi) 

•s!r"Vx)]l 

(43) 

In  the  above  equation,  the  Fresnel  reflection  coefficients  are  defined  as  follows: 

RTE{ku) 

_  fc\z  ^2  z 

k\z  "t  k2z 

(44) 

R™{ku) 

_  fyku  —  e\k2z 

(45) 

e2fciz  +  ex  k2z 

RTE(k2z) 

-  -RTE^  =  -t+t 

(46) 

R™(k2z) 

nTM ( /_  \  _  elklz  t\k2z 

—  K  [Kiz)  —  .  . 

^2  k\z  + 

(47) 

where  flTE(fciz)  and  R™(kiz)  are  the  reflection  coefficients  for  TE  and  TM  waves  inci¬ 
dent  from  region  Vi  to  region  V2,  respectively.  RTE(k2Z)  and  /?™(A;2z)  are  the  reflection 
coefficients  for  TE  and  TM  waves  incident  from  region  V2  to  region  Vi,  respectively. 


B.  Transmitted  Field  of  Equivalent  Source 

To  get  the  transmitted  wave  from  the  equivalent  source  in  the  lower  region  V2,  we  work 
on  the  up-going  wave  in  Eq.  (40)  as  following  steps: 

1.  Multiply  the  terms  with  e2  as  the  first  vector  in  the  tensors  with  Tte,  and  multiply  the 
terms  with  h2  as  the  first  vector  in  the  tensors  A2  and  B2  with  T™  , 
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2.  Change  the  first  vectors  e2(k2z)  and  h2(k2z)  in  the  polarization  tensors  to  ei(kiz)  and 
hi(kiz),  and 

3.  Change  the  up-going  wave  vector  k2  to  k\  in  the  exponential  terms. 

Thus  we  get 

m—  1 

/ dk±(ih,)m~'  (rrB(tj,)4i(*u)«j(*i.) 
+TTU(k2z)hl(kiz)h2(k2z))  •  a{rm\f±) 

+ 1 dk±(ik2,)m-'  (-T™(k2z)hl(ku)e2(k2,) 

+TTE(k2z)ei(ku)'h2(k2z))  •  6i“m)(fi)  J 

/ftiW1  (rTE(fc2Je,(/;u)e2(fc2l) 
+TTU(k2z)h,(kIZ)h2(k2,))  e^V-n 
•i  [V'x/(fl)  •  4"-m)(fl)] 

+  jdk±  (ik2,)m~2  (-T™ (k2z)hi(ku)e2(k2z) 


+TTE(k2z)et{ku)h2(k2z)] 

■i  [Vi/ft)  •  tf 

-"Vx)]} 

(48) 

where  the  transmission  coefficients  are  defined  as  follows: 

TTE(klz)  = 

'+RTE^=k?Xz 

(49) 

T™(klz)  = 

|(l  +  R™(kn)) 

772  2e2ku 

(50) 

f]i  e2kiz  +  €ik2z 

TTE(k2z)  = 

l+RTE^)  =  rrr 

Mz  “T  K2z 

(51) 

T™{k2z)  = 

l  (i + 
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=  &  *'k>\  ,  (52) 

r/2  e2fciz  +  eifc2z 

C.  The  First-Order  Field  Scattered  by  PEC  Rough  Surface 

For  a  perfectly  conducting  rough  surface,  the  tangential  electric  field  is  zero,  thus  = 
0.  Therefore  the  first-order  field  of  the  equivalent  source  can  be  written  as  [from  Eq.  (39)] 

J  dk±Ai  (kiz)e'*' (f _f,)  •  (fj.) 

K/(fl)  •  af(fl)]  .  (53) 

Notice  that  the  Fresnel  reflection  coefficients  of  the  TE  and  TM  waves  for  perfectly  con¬ 
ducting  surface  are  RTE(klz)  =  —1  and  R™(klz)  =  1.  Therefore,  from  Eq.  (43),  the 
first-order  (n  =  1)  reflected  wave  of  the  equivalent  source  from  the  perfectly  conducting 
surface  is 

^(f)  =  ^  j  J  dk±^ii(ku)ii(ku) 

+hi(klt)hl(ku)>j^<r-^  -  4Vx) 

+hl(klz)hl(kl2)y^<f-f,) 

^K/(fl)-ai0)(fl)] 

=  j df±f(r'±)  [ dk±~Ax{klz)jk'<f-fl)  •  af^) 

[Vl/(fl)-ai0)(fl)].  (54) 

Notice  that  the  first-order  reflected  wave  from  the  mean  surface  of  the  perfectly  conducting 
rough  surface  [Eq.  (54)]  is  the  same  as  the  first-order  radiated  field  of  the  equivalent  source 
[Eq.  (53)].  The  transmitted  field  is  zero,  thus  the  total  scattered  field  is 

eSV)  =  eSV) + ?(f)  =  2£j°(f) 
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=  —2 iki  J  dr'Lf(f,±)  ^  {Gt  (r,  r±)  •  a(°} (fj.) } 

+2^!  J  df'±Gi(f,f'x ) 

•*  K/(fl)  •  aiVx)] .  (55) 

IV.  Comparison  with  Conventional  SPM  Result 

Consider  a  plane  wave  Ei(r )  =  Ej°  (r)  =  ej£'0e**i‘f  incident  upon  a  perfectly  conducting 
rough  surface.  The  SPM  solution  of  the  zero-th  order  surface  variable  is  [7] 

ai0)(fl)  =  ai0)(fca)e<^^  (56) 


where 


d{x\ki±)  =qiaf)(ki±)  +piaf)(ki±) 
40)(M  =  2[ei-e(-kiz)]E01^ 

af\ki±)  =  2  [eV />(-*,*)] 


and  the  basis  vectors  are  defined  as  follows: 


Qi  =  —{xkiy  -  ykix) 
&ip 


Pi  —  Z{  x  ^  —  {xk{X  +  yk{y) 

&ip 

i{  =  5. 


Defining  the  Fourier  transform  pair 


and  the  delta  function 


F(ix)  =  (2^5 

/«.)  =  j  dk’^'^Ffc) 


S^  =  B wf*'-*1* 


we  can  write  the  following  relations  as: 

Wf±)  ,  r 


=  i  J  dk'±e«liF(Ic'±) 

=  i  j 

VJ(r i)  =  *  /  (k'xx  +  k'yy)  F(k'x). 
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By  making  "se  of  the  vector  products 


kx  ^iy  kyhix 


(xkx  +  yky)  •  pi  =  - 
( xkix  +  ykiy)  •  qi  =  0 
(xkix  +  ykiy)  -pi  =  k 

ei(kiz)-Qi  =  - 


kip 

kxkix 

+  ky^y 

kip 

=  0 

“  kip 

k*  k 
"’t  y^y 

+  kixkx 

kpk{p 

k^ky 

kiykx 

and  integrating  the  transverse  wavenumber,  the  first-order  scattered  field  in  Eq.  (55) 
rewritten  as 


kixky  k{ykx 


B'/’tf)  =  E°J dk±e‘l,  '<F(h  -  kii)ii(k  1.) 

\  rCpKip  #1  / 

+EC  j  dk±eil'-fiF(k±-ki±)ei(ku) 

+E„  j  dk^'HFik^-kiShxiku) 

(- 1 [a. .  e(— *,,)] 

\  kizkp  kip  J 

+E0 1 dkLerk'  fiF{k±  -  kix)hi(klz) 

(  O  ^  1  kxkix  +  kykiy  2  kpkip\ 

V  bp  klz  ) 

[e,-  •  /*(-fci2)]  . 

By  rewriting  the  scattered  electric  field  in  a  matrix  form 

-  [521 


L4l)(f)J 

J dh_jl'-fiF{kx  -  ki± ) 

I"  fee  feh  1  f  kiz)E0 


he  fhhl  LV  h{-kiz)E0 
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the  scattering  matrix  is  ‘hen  obtained  as 

fee  feh 

-  fhe  fhh  J 

■  qjl  kjvky-\-kjxkx  kjz  n l  kjxky  —kjvkx 

ZKl  kpkip  ki  ZKl  kpkip 

_ 2ki  kjz  kxkjy  kykjx  _  2k^  kx kjx~\-ky  kjy  j  2 kpkjp 

kitkp  k{p  k\tkp  k{p  k  1* 


(68) 


Expressing  the  wavenumber  components  by  using  angles  9,  <f>,  9,,  and  fa,  we  find  the  scat¬ 
tering  matrix  to  be 

fee  feh 
.  fhe  fhh . 

'  2 kiz  cos  (4>  —  <f>i)  2k\  sin  ( <j>  —  fa) 

.  ~%ki*  S^n  ~  <&)  ^  [Sin  ^  S^n  &i  ~  C0S  (0  ~  &)]  - 

(69) 

Notice  that  the  components  of  the  scattering  matrix  [Eq.  (68)]  derived  above  using  the 
equivalent  source  formulation  for  a  plane  incident  wave  are  the  same  as  the  ones  using 
conventional  SPM  as  given  in  [8]  for  the  limiting  case  of  k2  — >  oo. 


V.  Numerical  Result 


In  this  section,  the  backscattering  coefficients  will  be  calculated  by  using  the  first-order 
equivalent  source  formulation  for  a  Gaussian-tapered  incident  wave  and  compared  with 
the  standard  method  of  moments.  The  backscattering  coefficient  is  defined  as 


crab  =  lim  47rr2 
r— ►oo 


(l El  -  (El)  |2) 

Pi 


(70) 


where  Pi  —  f  |f?£|  dA  is  the  power  of  the  incident  wave  on  the  illuminated  area  with 

the  integral  over  the  mean  surface.  On  the  mean  surface,  the  electric  field  of  a  Gaussian- 

2/2 

tapered  wave  can  be  written  as  El  =  E0e~r  >s  ,  where  g  is  a  factor  specifying  the  beam 
width.  Therefore  we  get,  for  the  Gaussian-tapered  incident  wave, 


Pi  = 


00  2n 


J  rdr  J  #|£0|2e-2r’A2  =  \E0\ 


2  *9 


(71) 
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In  the  numerical  calculation,  the  expected  backscattering  coefficient  are  calculated  by 
averaging  the  samples  of  the  scattered  field,  thus 


crab  =  lim  A-kt4 

r-¥oo 


Pi 


47 rr2  1 

lim  — — -  — — - 
r-foo  p.  N  —1 

4?rr2  N 


N 


EW-j 


n=l 


21 


=  lim 


r-foo  P{  ft  _  1 


N 


iV  n=  1 


T.K 

n=l 


1  N 


(72) 


where  iV  is  the  number  of  realizations  in  the  Monte  Carlo  simulation  and  means 
statistically  estimated.  It  is  convenient  if  we  use.  the  standard  deviation  to  normalize  the 
calculated  profile  f(x,y).  The  standard  deviation  a  can  be  obtained  as 


=  P  kdk  daW(k,  a) 
J  o  Jo 


(73) 
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Fig.  4.  Comparison  of  the  scattered  field  calculated  by  using  the  spatial-domain  formulation  with  the 
result  by  the  method  of  moments  (MoM). 

In  the  numerical  simulation,  we  use  the  power  law  spectrum  to  generate  random  rough 
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surfaces.  The  power  law  spectrum  function  for  rough  surfaces  is  defined  as 

=  max 

*  (*max  -  kl)  k4 

where  a  is  the  rms  height,  and  kc  and  fcmax  are  the  low  and  high  cutoff  wavenumber  of  the 
rough  surface,  respectively. 

Fig.  4  is  the  numerical  result  for  the  backscattering  coefficients  of  the  co-polarized  waves 
<jhh  and  avv  by  using  the  spatial-domain  formulation  in  comparison  with  the  standard 
MoM.  The  rms  height  of  the  rough  surface  is  a  =  0.01A  and  the  number  of  realizations 
is  N  =  30.  The  incident  wave  is  tapered  with  the  factor  g  =  3.0A.  The  central  plane 
wave  of  the  tapered  incident  wave  varies  its  polar  angle  from  20°  in  order  to  avoid  the 
strong  specular  reflection  of  the  incident  beam  close  to  the  normal  incident  angle.  The 
comparison  of  the  results  demonstrates  good  agreement  between  the  two  different  methods. 
At  low  grazing  angles,  the  MoM  code  does  not  provide  correct  results  due  to  the  difficulty 
to  increase  the  size  of  rough  surface  at  these  angles  [9],  The  spatial-domain  formulation 
gives  reasonable  results  at  low  grazing  angles  where  ahh  drops  to  zero  and  crvv  tends  to  a 
constant  when  the  grazing  angle  approaches  zero. 


VI.  Discussions 

In  this  paper,  we  have  derived  the  spatial-domain  formulation  of  the  equivalent  sources 
for  any  orders  and  the  simple  spatial-domain  formulation  of  the  first-order  field  scattered 
by  a  PEC  rough  surface.  Up  to  the  first  order,  we  have  proved  that  the  total  field  of 
the  equivalent  source  on  a  PEC  rough  surface  is  the  field  scattered  by  the  rough  surface. 
The  advantage  of  using  the  spatial-domain  formulation  instead  of  frequency  domain  as 
in  the  conventional  SPM  formulation  is  that  there  is  no  need  of  integration  over  the 
spatial  frequency  and  the  non-averaged  complex  scattered  field  is  well-defined.  However 
the  derivation  for  the  total  scattered  field  of  higher  order  (n  >  2)  or  the  field  scattered 
from  a  dielectric  rough  surface  requires  much  more  effort  since  the  Sommerfeld  integrals 
are  involved.  By  considering  a  plane  incident  wave,  we  have  found  that  the  total  field  of 
the  equivalent  source  is  the  same  as  the  field  scattered  by  a  rough  surface  obtained  by 
using  the  conventional  small  perturbation  method  (SPM).  Particularly  for  the  PEC  rough 
surface,  the  equivalent  source  formulation  in  Eq.  (55)  is  more  suited  than  the  conventional 
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SPM  in  the  calculation  of  the  scattered  field,  since 

1.  The  incident  wave  E{  is  not  restricted  to  be  a  plane  wave  as  in  SPM,  therefore  it  works 
well  at  low-grazing  angles. 

2.  The  spatial-domain  formulation  can  be  used  to  calculate  the  near  field. 

3.  The  non-averaged  complex  scattered  field  in  spatial-domain  formulation  does  not  use 
the  Fourier  integral  over  the  rough  surface,  therefore  it  is  well-defined. 

4.  The  spatial-domain  formulation  is  computationally  efficient,  since  there  is  no  integration 
over  the  spatial  frequency  as  in  the  convention  SPM. 
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Appendix 

The  integral  representation  of  the  dyadic  Green’s  function  in  unbounded  space  is  given 
by  [7] 

+  — 2-  fjdkxdky— 

[Khz)e{ku)  +  h(ku)h(klz)]  eikl<f-^ 

Z  >  2? 

[e(-/=iz)e(-A;lz)  +  h(-ku)h(-kiz)] 

z  <  z* 

(75) 


where 


e(±klz) 

h{±ku) 

h 


xky  ykx 

Jk l  +  kl 


(xkx  +  yky) 


kxx  +  ky$  +  kltz 


(76) 


(77) 

(78) 
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K i  =  kxx  +  kyy  -  kizz 

(79) 

kn  =  sfk\  -  kl  - 

(80) 
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Electromagnetic  Wave  Interaction  of  Conducting  Object  with  Rough  Surface  by 

Hybrid  SPM/MoM  Technique 


I  Introduction 


Recently,  there  has  been  a  great  interest  of  studying  the  electromagnetic  wave  scattering 
from  object  situated  above  rough  surface  [1,  2,  3,  4,  5,  6,  7].  Previously,  electromagnetic 
wave  scattering  by  arbitrarily  shaped  objects  in  free  space  has  been  well  developed  using 
wire  [8,  9,  10]  and  surface-patch  model  [11,  12,  13,  14,  15,  16,  17,  18].  The  theory  and 
numerical  approaches  associated  with  object  near  flat  interface  of  layered  medium  have  also 
be  studied  extensively  by  many  researchers  [26,  27,  19,  20,  21,  22,  7,  31,  32,  25].  However, 
the  consideration  of  the  interface  to  be  rough  surface  is  a  new  challenge,  and  little  work 
has  been  reported.  In  theory,  the  standard  MoM  can  be  used  to  solved  unknowns  both  on 
object  and  rough  surface  [4,  28,  29].  However,  the  discretization  of  the  rough  surface  can 
significantly  increase  the  computational  resource  compared  to  calculating  the  scattering 
from  object  alone,  therefore  little  literature  has  been  reported  on  the  study  of  scattering 
for  full  scale  gemetry. 

In  this  paper,  we  present  a  hybrid  SPM/MoM  technique  to  calculate  the  EM  scattering 
from  3-D  conducting  object  above  rough  surface.  Under  the  hybrid  technique,  the  Green’s 
function  and  surface  variables  are  expanded  in  terms  of  the  surface  height  function  on 
the  mean  surface,  and  the  electric  integral  equations  based  on  the  extinction  theorem 
and  the  surface  boundary  conditions  are  decomposed  into  different  orders.  Each  order  of 
equations  represents  a  scattering  problem  with  the  same  geometry  and  different  equivalent 
incident  wave,  so  that  it  can  be  solved  efficiently  by  using  the  dyadic  Green’s  function  in 
layered  medium.  The  separation  of  solution  into  different  orders  also  helps  us  identify  and 
characterize  the  individual  interaction  terms  between  the  object  and  the  rough  surface. 

II  Configuration  and  Formulation 

Consider  an  electromagnetic  wave  Ei(f)  incident  upon  a  perfectly  conducting  object  with 
arbitrary  shape  S i  above  a  rough  surface  ST  as  shown  in  Fig.  1.  The  upper  space  Vi  is 
air  which  is  specified  by  t\  and  /zi,  the  lower  space  Vi  is  a  homogeneous  medium  and 
characterized  by  c 2  and  /Z2-  The  rough  surface  profile  is  defined  by  the  surface  height 
function  z  =  /(f'x)  with  mean  surface  S0  on  the  x-y  plane. 
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Figure  1:  Configuration  of  the  problem. 

A  Electric  Field  Integral  Equations 

On  the  conducting  surface,  the  tangential  electric  field  is  zero,  thus  we  can  write  the  integral 
equation  for  electric  field  as  following, 
for  f  €E  Si, 

(Ei(f)  +  fa  dS'  •  [n1(f/)  x  +  V  x  Vi{f,t)  •  jni(f')  x  ^(f')]} 

JJ^  dS'Gx(f,  f')  ■  x  n0(f)  =  0.  (1) 

By  Extinction  Theorem,  we  get, 
for  f  €  V 2, 

Ei(f)  +  fl  dS’  )  •  fni(f')  x  Hx(r  )]  +  V  x  Gx(r,f )  ■  [nj(f')  x  Ex{? )]} 

J  J  So 

+iW i  If  dS%(r,f')  •  Jtf)  =  0,  (2) 

and  for  r  E  V i , 

II  dS'  {*«W?2(f ,  f')  •  [n2(f')  x  772(f')]  +  V  x  S2(f ,  f')  •  [n2(f')  x  E2(fO] }  =  0,  (3) 


where  G\  and  G2  are  dyadic  Green’s  functions  in  unbounded  regions  with  (ej,//i)  and 
(e2,//2)  respectively.  The  vector  h\{f')  denotes  a  local  normal  pointing  from  rough  surface 
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Sr  to  upper  region  Vx.  The  vector  n2(f')  is  opposite  to  ni(f'),  thus  n2(r')  —  —hx{f').  E\{f') 
and  are  electric  and  magnetic  fields  on  the  rough  surface  in  region  Vx,  while  E2(f') 

and  H2{f')  are  surface  fields  on  rough  surface  in  region  V2.  Ji(r')  is  the  induced  surface 
current  on  the  object.  If  region  V2  is  dielectric,  the  tangential  fields  are  continuous,  thus 


_  _ 

ni(f')  x  H1(ff)  =  n^f')  x  H2(f')  = 


(4) 


_  „  Ht*  — 

n\{f)  x  Ei(f')  =  ni(f')  x  E2{f')  = 


(5) 


where  a(rx)  and  b(r'x)  are  new  surface  variables  defined  on  the  mean  surface  S0,  ifo  is  the 
characteristic  impedance  of  upper  region  and  can  be  written  as  rft  =  (fXi/e i)1/2,  and 


dr'x 


dS‘  yjl  +  (df(x\y^)/dx^f  + 

With  the  new  surface  variables,  we  can  rewrite  the  integral  equations  as 

(Ei(r)  +  df'jt  {iAqGi  (f,  f')  •  d(r'x)  +  V  x  Gx  (f,  f')  ■  6(r^)} 
+iojpii  JJ  dS'Gi(f,r')  ■  Ji(f r)j  x  n0(f)  =0  for  f  G  Si 


(6) 


(7) 


Ei(r)  +  ff  {ifciGi  (r,r') -a(rx)  +  V  x  Gi  (r,r')  •  6(r^)| 

J  J  So 

+iu>ni  JJ  dS'Gi(f,r')-Ji(f')=  0  for  r  E  V2, 


1 f  dr'xUk2— 
JJso  {  r)  1 


G2  (r,  r')  ■  a( f^)  +  V  x  G2(r,  r')  •  b(rx)  }  =  0 


for  r  €  V\. 


(8) 

(9) 


Theoretical,  given  the  rough  surface  ST  and  the  perfectly  conducting  object  surface 
profiles  Si,  the  unknown  surface  variables  a(f'x),  b(fx)  and  the  induced  current  Ji(f')  can 
be  solved  from  Eqs.  (7-9).  For  the  special  case  in  which  the  interface  is  flat,  the  surface 
variables  d(rx)  and  b(fx)  are  in  the  horizontal  plane  and  the  local  coordinate  f'  can  be 
replaced  by  r'x  in  the  dyadic  Green’s  functions  Gi(f,  f)  and  G2(f,  f'). 
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B  Expansion  of  Green’s  Function  and  Surface  Variables 

We  assume  that  the  rough  surface  height  is  small.  Therefore  the  scalar  Green’s  function 
in  region  a  (where  a  =  1,  2)  can  be  expanded  in  terms  of  surface  height  function  /(f'x)  on 
the  mean  surface  SD, 


9a(f,f')  =  ga{f,f'±  +  zf(f'x)) 

co  i  am  co  i  om 


„0m!dzm' 


771=0 


m\  dzm' 


in  which  the  following  property  of  the  Green’s  function  has  been  used 

d  d  , 

Thus  the  dyadic  Green’s  function  can  be  written  as 


(ii) 


G«(f,  f)  =  (/  +  iw)  g„(f,  f)  =  £  T  (-/(rl))”1  (f,  fi) .  (12) 


Consequently,  the  surface  variables  and  induced  current  can  be  written  as  the  series 
expansions  as  following 

Ji(7)=  ETfV).  (13) 

771=0 

OO 

a(f±)  =  y,  [^rVi) + ^im)(^L)] .  (i4) 

777=0 

°° 

Krl)  =  E  RV±)  +  ^’(fl)]-  (IS) 

771=0 

In  Eqs.  (14)  and  (15),  the  separation  of  the  ^-components  for  the  surface  variables  a  and 
b  allows  the  unknown  surface  variables  become  only  the  tangential  components  since  the 
z-components  of  the  m-th  order  can  be  expressed  in  terms  of  the  tangential  components 
of  lower  order  (m— 1)  as  we  will  see  below.  By  the  definition  of  a  and  b  as  in  Eqs.  (4)  and 
(5),  we  get  the  following  identities, 


where 


nl(r1)-a(r,x)  =0, 

(16) 

ni{r±)  -  b(f’x)  =  0, 

(17) 

,  ,  V'±/(f'J  +  z 

\V’J(r'±)  +  z\ 

(18) 
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Substituting  the  series  expressions  of  a  and  b  as  in  Eqs.  (14)  and  (15)  into  Eqs.  (16)  and 
(17)  respectively,  and  considering  that  df{f'L)/dx'  and  df(r'±)/dy'  are  in  the  same  order 
as  we  then  get  the  m-th  order  z-components  of  surface  variables  a  and  b  in  terms 

of  the  (to  —  l)-th  order  tangential  components, 


=  V‘J( fL)  ■  af-'V-l), 


(19) 


Substituting  the  series  expression  of  dyadic  Green’s  functions  as  in  Eq.  (12)  and  the 
surface  variables  a  and  b  as  in  Eqs.  (14)  and  (15)  into  the  integral  equations  (7-9),  we  get 
the  n-th  order  equations  in  terms  of  the  surface  height  function  f(f'±)  as  following. 

C  The  n-th  Order  Equations 


d“V)  +  Jh  {ikiGi  (f,  f±)  ■  +  VXC,  (f,f i) :  5?>(fi)} 

+iu/ii  JJs  x  ho(f)  =  0  for  f  €  Si,  (20) 


B?\f)  +  JJ  d7±{iktGl(fyi)-di?\r'J  +  VxG^f,fll)-b™(f'±)} 


(21) 


-kiujfii  jJs  dS'Gi (r,  r')  -  J'i'iV)  =0  for  r  e  V2l 

S"’(f)+//s for  f  e  Vi, 
where  for  n  =  0, 


(22) 


£f(f)  =  E,(  f), 

Eg’(f)  =  0, 


(23) 

(24) 


for  n  >  1, 


dz' 


{ik,G,  (f.fV)  ■  af-Vj  +  V  x  G,  (f.fp  ■  S,;-"”(f'i)} 
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I’kiG,  (f, fi)  •  z  [vjif’,)  •  a("“'n,(f'1)] 

+Vxff1(f,f'1)l[v'1/(f'i)B'Ln-’n,(fi)]},  (25) 


dz” 


ik2 -G,  (f,  fi)  •  a(rm>(fi)  +  V  x  G2(f,  fi)  •  5<”-”>(fi) 


Vi 


jll/dfil— /(fi)]’ 


.-i  d 


m— 1 


dz 


m-1 


{ H=2^p2  (f.fi)  •  i  [Vi/(fi)  •  ai”-m)(fi)] 
+v  x  52(f,fi)  •  %  [vi/(fi)  •  5!rm)(fi)] } . 


(26) 


Comparing  the  zeroth  order  (n  =  0)  integral  equations  with  Eqs.  (7-9)  for  flat  interface 
(/(f±)  =  0)»  we  found  that  they  are  the  same  except  adding  a  superscript  (0)  to  the  surface 
variables  ai  and  f>j_  and  the  induced  current  J\ .  Therefore,  the  solutions  and 

for  the  zeroth  order  equations  should  be  the  same  as  the  ones  for  Eqs.  (7-9)  in  which  the 
rough  surface  is  considered  to  be  flat.  For  the  higher  order  equations  (n  >  1),  we  found  that 
they  are  also  in  the  same  form  as  the  zeroth  order  equations,  except  for  the  substitution  of 
Ei  by  E[n\  and  the  additional  “source”  term  E^  in  the  lower'region.  Therefore  we  only 
need  to  know  how  to  solve  the  zeroth  order  equations,  i.e.  the  equations  for  object  on  a 
flat  interface. 

Since  the  equations  of  any  order  are  equivalent  the  ones  for  flat  interface,  they  can  be 
rewritten  by  introducing  the  dyadic  Green’s  function  for  layered  media.  The  advantage  for 
this  approach  is  avoiding  solving  surface  unknowns  a ^  and  on  interface,  but  only 
need  to  solved  the  induced  current  on  the  conducting  body.  Therefore  both  the 
computational  time  and  memory  can  be  dramatically  saved. 
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D  Application  to  N on-penetrable  Rough  Surface 

Without  losing  the  generality,  we  consider  the  rough  surface  to  be  non-penetrable.  In  this 
case  the  zeroth  order  electric  field  integral  equation  (EFIE)  can  be  written  as 

h0(f)  x  ^E-0)(f)  +  E^ir)  +  iwfii  dS'GL(r,  f')  •  jf\r)J  =  0 

and  the  n-th  order  (n  >  1)  EFIE  can  be  written  as 

n„(f)  x  i  ^  dS'Gtif,  r)  •  jjn)(f/) j  =  0 

where 

=  E  (-/(fir dm 

m=l 


dzr 


for  r  €  Si, 

(27) 

for  f  €  Si, 

(28) 

n  1  r  r  *  ~  1 


dz™-1 

{ifcA  (f,  fl)  •  z  [Vx/(fj_)  •  4n-m)(fl)] }  , 
and  the  dyadic  Green’s  function  for  conducting  interface  is 

VL(f,  f')  =  Si  (f ,  f')  -  Si(f,  (7  -  2zz)  •  r')  ■  (7  -  2zz), 


(29) 


(30) 


where  G\  is  the  dyadic  Green’s  function  in  unbounded  medium  of  region  V\ . 

The  surface  field  dj_  of  lower  order  can  be  obtained  from  the  lower  order  scattered  field  from 
object  and  equivalent  source.  For  example,  the  first  order  surface  variable  a ^  can  be  calculated 
by  using  the  solution  of  the  zeroth  order  induced  current  and  the  incident  and  reflected 
magnetic  fields  in  absence  of  the  object,  T  .  ..  . 

4V l)  =  771  z  x  (Hn(f'±)  +  Hri(f'±)  +  jf  dS*Y  x  SL(f(L,f")  •  jfV))  •  (31) 

The  higher  order  (71  >  2)  surface  variable  0-4  can  be  obtained  similarly  but  involving  much  more 
manipulation. 

Up  to  the  first  order,  the  total  returned  field  can  be  written  as 

Es(f)  =  Er(f)  +  Eb(f)  +  Ec(f)  +  Ed(f),  (32) 

where 

Eb(f) 

Ec(f) 

Ed(f) 


=  iuifii  I  dS'Gi{r,r')  ■  (33) 

Js1 

=  E\l)(f)  +  E{rl\f),  (34) 

=  ium  [  dS,GL(f,r,)-j\l)(r').  (35) 
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Figure  2:  Radiation  and  reflection  of  the  equivalent  source. 


and 

£iX)(f)  =  ~^i  JJ  df'±f(f,x)-^-Gi  (f.f'J  •of)(f/x) 

+ifci  / f  df'jG i  (f,fl)  •  z  [v'x/(fi)  •  df  (#=!)] .  (36) 

In  Eq.  (32),  the  £r(f)  is  simply  the  reflected  field  from  the  flat  interface  in  absence  of  con- 
ducting  object.  In  the  expression  of  Eb(f)  as  in  Eq.  (33),  the  induced  current  J\  1  is  obtained 
by  solving  the  integral  equation  (27)  with  layered  Green’s  function,  therefore  the  returned  field 
Efc(f)  includes  all  interactions  between  object  and  flat  interface.  The  Ec{f)  as  in  Eq.  (34)  is  the 
sum  of  radiation  field  from  the  “equaivalent  source”  and  its  reflection  are  illustrated  in  Fig.  2. 
The  reflected  field  of  the  equivalent  source  can  be  obtained  by  writing  unbounded  dyadic  Green’s 
function  in  the  integral  form  as  in  Appendix  A,  thus  the  radiation  field  of  the  equivalent,  source 
e[  ‘  is  repressed  as  the  sum  of  plane  waves.  By  knowing  the  Fresnel  reflection  coefficient  of 
each  plane  wave  component,  it  is  easy  to  write  the  reflected  field  E^  by  multiplying  RTE  and 
R™  to  the  reflected  TE  and  TM  waves  respectively.  For  non-penetrable  surface  RTE  =  —  1  and 
rTM  _  ^  thus  we  find 

4l)(f)=^(1)(f).  (37) 

Therefore  the  return  field  due  to  the  equivalent  source  is  simply 

Ec(f)  =  2E(il\f).  (38) 

It  can  also  be  shown  that  the  returned  field  Ec(f)  as  in  Eq.  (38)  is  the  same  as  the  first 
order  SPM  solution  for  conducting  rough  surface  if  we  let  the  induced  current  =  0  when 
evaluating  the  surface  field  a^(f^)  in  Eq.  (31).  Therefore,  we  call  the  return  field  Ec(r)  to  be 
the  ’’incoherent”  returned  field  from  rough  surface  under  the  influence  of  the  object.  The  returned 
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field  £</(f)  is  the  radiation  field  of  the  first  order  induced  current  excited  by  the  ’’incoherent” 
field  f?c(r).  Since  the  layered  Green’s  function  are  used  to  calculate  the  induced  current  j[1^  and 
its  radiation  field,  therefore  the  returned  field  Ed{f)  includes  all  multiple  interactions  between 
object  and  conducting  interface. 


Ill  Numerical  Results 


In  the  numerical  simulation,  a  horizontal  conducting  cylinder  along  x-axis  with  2.0A  in  length  and 
1.0 A  in  diameter  is  considered.  The  distance  between  the  bottom  of  the  cylinder  and  the  mean 
height  of  conducting  rough  surface  is  0.1  A  so  that  large  interaction  between  object  and  rough 
surface  can  be  involved.  A  rough  surface  with  well-known  Gaussian  power  spectrum  is  used  for 
the  validation  with  the  standard  method  of  moment  (MoM).  The  size  of  the  rough  surface  is  15. 0A 
by  15.0A.  The  deviation  and  correlation  length  of  the  rough  surface  are  a  =  0.03A  and  l  =  1.0A 
respectively.  The  incident  wave  E{  is  tapered  and  formed  as  the  summation  of  plane  waves  with 
Gaussian-shaped  footprint  on  the  mean  surface.  The  tapered  incident  wave  satisfies  Maxwell’s 
equations  and  minimizes  the  edge  effect  in  numerical  calculations.  The  factor  g ,  which  is  used  to 
control  the  beam  width  of  the  tapered  wave  is  g  =  3.0A,  so  that  the  incident  electric  field  density 
on  the  rough  surface  drops  1/e  at  \rfx\  =  3.0A. 

In  the  numerical  calculation,  the  radar  cross  section  (RCS)  is  defined  as 


RCS  =  lim  Anr2 


IMM)f 


(39) 


\E0{Qi,4>i)\2' 

where  E0  is  the  maximum  magnitude  of  tapered  incident  electric  field  on  the  mean  surface  SQ ■ 
For  monostatic  (backscattering)  RCS,  the  scattering  angles  6  ~  6{  and  <p  =  <f>{.  In  the  numerical 
simulation  for  bistatic  RCS,  we  let  the  scattering  6  angle  be  40°  and  vary  the  azimuthal  angle  <p 
from  0°  to  360°,  hence  <j>  =  0°  is  the  backscattering  direction.  For  monostatic  RCS  simulation, 
the  incident  k{  vector  remains  in  the  x  —  z  plane  and  the  incident  angle  0{  varies  from  0°  to  90°. 

The  plots  shown  in  Fig.  3  and  Fig.  4  are  bistatic  RCSs  for  individual  terms  of  Eq.  (32)  for  TE 
and  TM  incident  waves  respectively.  The  plot  labeled  as  Er  is  the  reflected  field  of  the  incident 
tapered  wave  from  flat  interface  in  absence  of  object.  A  peak  of  co-pol  component  of  Er  appears 
in  the  specular  direction  <f>  =  180°  as  expected.  The  plot  labeled  as  E{,  is  the  returned  field  as  in 
Eq.  (33).  The  Ec  plot  is  evaluated  by  using  Eq.  (34),  and  the  Ed  plot  is  calculated  from  Eq.  (35). 
It  can  be  found  that  most  of  the  energy  of  the  “incoherent”  field  Ec  concentrates  in  the  forward 
scattering  direction.  The  fields  shown  in  the  E &  plot  are  the  scattering  from  object  excited  by 
fields  Et  and  Er ,  while  the  secondary  scattering  field  in  the  Ed  plot  is  the  returned  field  from  the 
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object  excited  by  the  “incoherent”  field  Ec  from  the  rough  surface.  We  note  that  the  cross-pol 
returned  fields  VH  and  HV  in  the  backscattering  direction  are  zero  in  the  E &  plot.  This  is  due  to 
the  symmetrical  property  of  the  geometry  and  the  incident  wave.  The  secondary  returned  fields 
of  cross-pol  VH  and  HV  in  the  plot  E 4  are  no  longer  zero  in  the  backscattering  direction  due  to 
the  asymmetry  of  the  “incoherent”  field  Ec.  The  “incoherent”  field  Ec  as  well  as  the  secondary 
returned  field  Ed  from  the  object  are  both  proportional  to  the  surface  height  function  of 

the  rough  surface.  It  can  be  easily  checked  that  both  Ec  and  Ed  become  zero  when  the  surface 
height  function  is  zero. 

The  sum  of  the  four  terms  Er  +  Eb  +  Ec  +  Ed  as  in  Eq.  (32)  is  the  total  returned  field  up  to  the 
first  order  as  shown  in  Fig.  5  and  Fig.  6  for  TE  and  TM  incident  wave  respectively.  The  bistatic 
RCS  for  the  total  returned  field  with  TE  incident  wave  are  compared  with  standard  MoM  results. 
In  the  MoM  simulation,  both  the  rough  surface  and  the  conducting  object  are  discretized  and  the 
surface  unknowns  are  solved  together  by  using  the  conjugate  gradient  algorithm.  The  plots  in 
the  top  columns  in  Fig.  5  and  Fig.  6  show  the  bistatic  RCS  of  the  zeroth  order  solution  using  the 
hybrid  technique  and  the  standard  MoM  solution.  Good  agreement  is  found  in  the  comparison  of 
the  results.  The  co-pol  W  of  the  zeroth  order  solution  matches  better  than  the  co-pol  HH,  since 
the  TE  incident  wave  may  induce  unnegligible  currents  on  the  front  and  back  edges  of  the  rough 
surface,  which  may  add  more  returns  for  the  MoM  result.  This  edge  effect  can  be  minimized  if 
the  TM  wave  is  used  as  incident  as  we  can  see  in  Fig.6.  In  the  simulations  as  shown  in  Fig.  5 
and  6,  a  single  rough  surface  with  Gaussian  power  spectrum  is  used.  The  deviation  of  the  surface 
height  is  cr  =  0.03 A  and  the  correlation  length  is  l  =  1.0A.  The  results  also  show  good  agreement 
in  comparison  of  the  results  obtained  by  hybrid  technique  with  the  standard  MoM.  It  is  noted 
that  the  curves  are  no  long  symmetrical  with  respect  to  the  incident  plane  when  the  surface  is 
rough. 

The  Monte  Carlo  simulation  results  with  100  realizations  are  shown  in  Fig.  7  and  Fig.  8  for  TE 
and  TM  incident  waves  respectively.  The  rough  surfaces  with  power  law  spectrum  are  generated 
and  they  are  independent  to  each  other.  The  power  law  spectrum  function  is  W(k)  ~  a0/kA  which 
more  closely  represents  ocean  state,  where  k  is  the  spatial  wavenumber  of  the  rough  surface  and 
a0  =  0.008/27T  which  corresponds  to  the  amplitude  used  for  Durden-Vesecky  spectrum  [35].  The 
higher  cut-off  spatial  wavenumber  k ^  is  chosen  to  be  k h  =  2.5/ci  which  corresponds  the  band  width 
for  l/5Ai  spatial  resolution  of  sampling  on  the  rough  surface,  where  k\  and  Ai  are  electromagnetic 
wavenumber  and  wavelength  of  incident  wave  respectively.  The  lower  cut-off  spatial  wavenumber 
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MoM 

Hybrid 

Number  of  unknowns 

17,252 

972 

CPU  time 

single  bistatic 

100  realizations 

single  monostatic 

100  monostatic 

1  hr  14  min 

5  days  (estimated) 

23  hr 

3  months  (estimated) 

5  min 

40  min 

3  hrs  15. min 

17  hrs 

Table  1:  Comparison  of  computational  effort  for  standard  MoM  and  hybrid  codes,  respec¬ 
tively. 

ki  is  chosen  according  to  the  deviation  of  rough  surface  height  using  the  following  relation 

<40) 

In  Fig.  7  and  Fig.  8,  the  deviations  of  rough  surface  k\a  =  0.1,  0.2,  0.4  correspond  to  the  lower 
cut-off  spatial  wavenumbers  k{  ~  0.6131,  0.3137,  0.1578fci  respectively.  It  is  shown  that  the 
Monte  Carlo  simulation  converges  with  respect  to  the  number  of  realizations,  and  the  averaged 
cross-pol  RCSs  increase  with  the  deviation  of  rough  surface  at  a  wide  range  of  scattering  angles. 

Fig.  9  and  10  show  the  Monte  Carlo  simulation  of  monostatic  (backscattering)  RCS  for  TE  and 
TM  incident  waves  respectively.  Here  100  rough  surfaces  with  power  law  spectrum  and  k\a  =  0.4 
are  used.  The  backscattering  direction  varies  from  6  =  0°  to  6  =  90°  with  45  steps  in  between. 
The  azimuthal  angle  <p  remains  0  degree.  We  note  that  the  cross-pol  VH  and  HV  are  significant  in 
the  presence  of  rough  surface.  Analytically  the  cross-pols  should  be  zero  for  the  geometry  we  are 
considering  when  the  surface  is  flat.  We  note  that  the  non-zero  cross-pols  for  flat  surface  shown 
in  the  solid  lines  are  produced  numerically.  The  changes  of  co-pol  monostatic  RCS  due  to  the 
rough  surface  can  be  found  at  some  scattering  angles.  At  small  grazing  angle,  the  rough  surface 
effect  on  co-pol  backscattering  RCS  is  not  significant. 

The  computational  performance  of  the  hybrid  technique  and  the  standard  MoM  are  listed 
in  Table  1  for  the  testing  cases  discussed  above,  in  which  the  Alpha  21164  machine  with  2MB 
L2  Cache,  500  MHz  clock  rate  and  1GB  RAM  are  used  for  the  simulation.  The  CPU  time  for 
monostatic:  simulations  are  based  on  45  incident  angles. 


=  J  W(k)d2k  = 


7T  CLq 
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IV  Conclusion 


This  paper  presents  a  hybrid  technique  of  SPM  and  MoM  for  EM  scattering  from  object  above 
rough  surface.  With  the  expansion  of  the  Green’s  function  and  surface  variables  in  terms  of  the 
surface  height  function  on  the  mean  surface,  the  electric  field  integral  equations  are  decomposed 
into  different  orders.  The  equations  of  each  order  represent  the  EM  scattering  problem  with  the 
same  object  above  the  mean  surface  (flat)  and  different  incident  field  from  equivalent  source  which 
can  be  evaluated  by  using  lower  order  solutions.  The  equivalence  of  the  flat  surface  problem  allows 
us  to  use  the  dyadic  Green’s  function  for  layered  media,  so  that  we  do  not  need  to  solve  for  tan¬ 
gential  fields  on  the  rough  surface,  leaving  only  unknowns  on  the  conducting  object.  Comparison 
with  standard  MoM,  this  hybrid  technique  demonstrates  a  dramatic  increase  in  computational 
efficiency  without  losing  the  accuracy.  The  separation  of  return  field  into  the  sum  of  individual 
interaction  terms  allows  us  to  identify  the  coherent  and  incoherent  returned  field,  therefore  the 
rough  surface  affect  can  be  quantitatively  characterized. 
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Appendix  A 


The  integral  representation  of  dyadic  Green’s  function  in  unbounded  space  is 


G\{r,r  )  =  -zz 


+ 


5 

^  ff  dkzdky^ 


e(klz)e(klz)  +  h(klz)h(klz)}  e'Mf-f' 
e(-ku)e(-kiz)  +  e 


z  >  z\ 
iKV(f-f')  z<z!y 


(41) 


where 


e(±kiz) 

h(±klz) 


xky  ykx 

AX  +  k^ 


'-lz 


yfk. 


I  + 


(xkx  +  yky )  +  z 


(42) 

(43) 


k\  —  kxx  kyfj  “f-  k\zz 


(44) 


13 


hxi  "f-  y  1 2 2 


(45) 

(46) 


K  i  = 

k\z  —  \j^{  —  k^.  —  ky 
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Figure  4:  Bistatic  RCS  of  individual  terms  for  TM  incident  wave. 
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Figure  6:  Bistatic  RCS  of  the  total  returned  field  Er  +  Eb  +  Ec  +  Ed  for  TM  incident  wave. 
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Figure  7:  Monte  Carlo  simulation  of  bistatic  RCS  of  the  total  returned  field  Er+Eb+Ec+Ed 
for  TE  incident  wave. 
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Figure  8:  Monte  Carlo  simulation  of  bistatic  RCS  of  the  total  returned  field  ET+Eb+Ec+Ed 
for  TM  incident  wave. 
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Figure  9:  Monte  Carlo  simulation  of  monostatic  RCS  of  the  total  returned  field  Er  +  Eb  + 
Ec  +  Ed  for  TE  incident  wave. 


Figure  10:  Monte  Carlo  simulation  of  monostatic  RCS  of  the  total  returned  field  Er  +  Eb  + 
Ec  +  Ed  for  TM  incident  wave. 


23 


Tapered  Wave  with  Dominant  Polarization  State  for  All  Angles  of  Incidence 


2 


Abstract 

Typical  applications  of  the  method  of  moments  (MoM)  to  rough  surface  3-D  electromagnetic  scattering 
require  a  truncation  of  the  surface  considered  and  call  for  a  tapered  incident  wave.  It  is  shown  how  such 
wave  can  be  constructed  as  a  superposition  of  plane  waves,  avoiding  problems  near  both  normal  and 
grazing  incidence  and  providing  clean  footprints  and  clear  polarization  at  all  angles  of  incidence.  The 
proposed  special  choice  of  polarization  vectors  removes  an  irregularity  at  the  origin  of  the  wavenumber 
space  and  leads  to  a  least  squared  error  property  of  the  wave.  Issues  in  the  application  to  3-D  scattering 
from  an  object  over  a  rough  surface  are  discussed.  Approximate  3-D  scalar  and  vector  tapered  waves  which 
can  be  evaluated  without  resorting  to  any  numerical  integrations  are  derived  and  important  limitations 
to  the  accuracy  and  applicability  of  these  approximations  are  pointed  out. 

Keywords 

Tapered  wave,  rough  surface  scattering,  method  of  moments. 

I.  Introduction 

Recent  years  have  seen  major  advances  in  the  development  of  fast  MoM  solvers  for  3-D 
scattering  of  electromagnetic  vector  waves  from  rough  surfaces  [1],  [2],  [3],  [4],  [5],  [6], 
[7],  [8].  Efforts  are  now  also  being  directed  towards  inclusion  of  objects  situated  in  the 
neighborhood  of  the  rough  surface  [9],  [10],  [11],  [12].  Since  the  problem  of  scattering 
from  an  object  next  to  a  rough  surface  is  computationally  complex,  2-D  investigations  are 
also  of  importance  [13],  [14],  [15],  [16],  [17],  [18],  [19],  [20],  The  3-D  case  with  or  without 
objects  is  aimed  at  by  the  present  paper.  •'  - 

The  methods  employed  usually  require  a  truncation  of  the  rough  surface  because  of  lim¬ 
ited  computing  resources  which  leads  to  erroneous  results  due  to  artificial  edge  diffraction 
when  ideal  plane  waves  are  used  to  excite  the  system.  The  tapered  wave  concept  is  based 
on  providing  an  illumination  for  the  numerical  simulation  that  resembles  the  plane  wave 
case  to  be  modeled  closely  at  the  center  of  the  scattering  scenario  (including  a  particu¬ 
lar  arbitrary  polarization)  while  its  intensity  becomes  negligibly  small  upon  approaching 
the  artificially  introduced  edges  of  the  rough  surface.  Thus  unwanted  edge  effects  due  to 
the  primary  incident  wave  are  avoided  and  the  proper  normalization  of  computed  scatter¬ 
ing  coefficients  allows  a  meaningful  comparison  with  the  ideal  plane  wave  case;  near-field 
quantities  such  as  current  distributions  induced  near  the  center  of  the  tapered  wave  are 
also  expected  to  be  similar. 
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Furthermore,  the  tapered  wave  should  be  constructed  in  such  a  way  that  it  satisfies  the 
Maxwell  equations  without  any  approximation.  This  helps  to  increase  the  confidence  in 
the  results  obtained  from  the  in  general  rather  complex  MoM  simulation  codes.  It  should 
also  be  possible  to  substitute  it  for  a  plane  wave  of  arbitrary  polar  and  azimuthal  angles 
of  incidence  without  loss  of  polarization  and  degradation  of  tapering. 

The  above  requirements  led  us  to  revise  and  modify  the  tapered  wave  found  in  the  open 
literature  which  is  based  on  a  superposition  of  plane  waves. 

II.  Superposition  of  Plane  Waves 

Consider  a  homogeneous,  isotropic  medium  with  real  wavenumber  k  and  wave  impedance 
77.  Then  the  following  superposition  of  a  2-D  spectrum  of  plane  waves  is  an  exact  solution 
to  the  Maxwell  equations  and  represents  a  wave  incident  upon  the  x-y  plane  from  z  >  0: 

OO 

Ei(f)=  J  dkpe^-^z)  ^(kp)e(kp)  (1) 

—  OO 

Hi(r)  =  7  dkpei{lp  ~p-k^z)^^-h(kp)  (2) 

J  Tj 

—  OO 

where 

r=p  +  zz  (3) 

kp  =  x  kx  +  y  ky  .... .  .  ...  (4) 

and 


The  spectrum  ip  carries  the  information  on  the  shape  of  the  footprint  (defined  as  the 
distribution  of  the  magnitude  in  the  x-y  plane)  of  the  incident  field  and  also  on  the 
direction  of  incidence.  It  is  assumed  to  be  centered  about 

k{p  x  kix  -P  y  kiy  (b) 

=  k  sin  6i  (x  cos  (pi  +  y  sin  (pi)  (7) 

where  6i  and  (pi  are  the  polar  and  azimuthal  angles  of  incidence  of  the  central  plane  wave 
and — pars  pro  toto — of  the  tapered  wave.  In  an  application  the  central  plane  wave  would 


kz  kz(kp)  \ 


J* 


k2  0  <kp<  k 


A  p  ^  A/ 


(5) 
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coincide  with  the  pla*'  wave  which  was  replaced  by  the  tapered  wave  in  the  numerical 
simulation.  Details  about  the  functional  dependence  of  ip  are  given  in  Section  III. 

The  polarization  vectors  e  and  h  are  of  the  general  form 

e(kp)  =  eh(kp)  h(kp)  +  ev(kp )  v(kp )  (8) 


h(/up)  —  o/i(&p)  f^(^p) 


The  notations 


h{kp)  = 


x  sin  &  —  y  cos  4>i  kp  =  0 


,  ( x  ky  y  k% ) 


kp  >  0 


v(kp)  = 


x  cos  4>i  +  y  sin  (pi 


kp  —  0 
kp  >  0 


[  +  2/M  +  ^  if  >  0 

are  found  in  similar  form  in  [21],  [22].  The  chosen  definitions  for  kp  =  0  take  care  of  the 
special  case  of  a  normally  incident  pure  plane  wave.  kp  >  k  corresponds  to  evanescent 
waves  (Section  III)  and  the  horizontal  part  of  v(kp)  is  imaginary  in  this  case.  It  is  important 
to  note  the  discontinuity  of  h  and  v  at  kp  =  0;  both  unit  vectors  change  sign  when  crossing 
the  origin  along  a  straight  line  in  the  kp  plane. 

The  general  superposition  integrals  (1),  (2)  were  stated  similarly  in  [23],  [7];  however, 
only  normal  incidence  is  considered  in  what  follows  there.  If  e/,  and  ev  in  (8),  (9)  are  set 
to  constants  then  (1),  (2)  specializes  to  the  tapered  wave  used  in  [1],  [2],  [4],  [9],  [8],  with 
a  particular  spectrum  ip  briefly  discussed  in  Section  III.  Problems  with  this  tapered  wave 
encountered  near  the  grazing  incidence  (for  discussion  and  references  see  Section  III)  and 
near  the  normal  incidence  (Section  IV)  motivated  our  formulation  of  a  different  kind  of 
tapered  wave,  especially  with  respect  to  the  polarization  vectors. 

III.  Amplitude  Spectrum 

If  the  polarization  vector  e  on  the  right-hand  side  of  (1)  is  replaced  by  a  scalar  constant 
then  the  resulting  integral 


£i(r)  =  e,  J  dkpe 1 


(kp'P—kzz) 


ip(kp) 


April  24,  2000 


DRAFT 


5 


is  the  plane-wave  representation  of  a  scalar  wave  satisfying  the  scalar  Helmholtz  equation 
and  ip  can  be  identified  with  the  well-known  angular  spectrum  in  scalar  diffraction  theory 
[24],  [25].  Thus  by  obtaining  ip  via  2-D  Fourier  transformation  and  making  sure  that  e 
and  h  vary  only  moderately  over  the  spatial  frequency  range  where  ip  is  not  negligible, 
arbitrary  footprints  of  the  vector  tapered  wave  can  be  approximated.  [The  mentioned 
requirement  leads  to  a  problem  with  the  tapering  in  [1],  [2],  [4],  [9],  [8]  near  normal 
incidence  (Section  IV).]  The  information  on  the  direction  of  incidence  of  the  tapered  wave 
is  included  by  shifting  ip  in  the  kx-ky  plane  to  be  centered  about  kip.  The  prescribed 
footprint  itself  is  fixed  with  respect  to  angle  of  incidence. 

A  Gaussian-shaped  footprint  (Fig.  1)  whose  amplitude  at  p  =  g  is  down  to  1/e  times 
the  level  at  the  center  is  implemented  by  choosing 

i>(kp)  =  e (13) 

A  pure  plane  wave  is  described  by  ip(kp)  =  S(kp  -  kip),  or  g  -»  oo  in  (13).  It  should 
be  pointed  out  that — as  is  well  known  from  signal  theory — among  all  footprints  of  given 
finite  energy  and  width,  the  Gaussian  leads  to  the  smallest  bandwidth  (for  the  appropriate 
definition  of  space  and  frequency  domain  widths)  which  is  desirable  for  synthesis. 

Spectral  components  with  kp  >  k  are  the  amplitudes  of  plane  waves  that  travel  along 
the  x-y  plane  and  are  evanescent  for  z  >  0.  Their  inclusion  makes  it  possible  to  synthesize 
a  given  footprint  near  or  at  grazing  incidence. 

-  The  spectrum  in  [1],  [2],  [4],  [9],  [8]  is  given  as  a  2-D  Fourier  integral  that  needs  to  be 
evaluated  numerically.  It  cannot  be  used  near  grazing  incidence  where  the  field  distribution 
in  the  x-y  plane  becomes  highly  oscillatory.  Its  continued  use  is  rooted  in  its  close  relation 
to  a  scalar  tapered  wave  employed  previously  [26],  [27].  The  latter  wave,  on  the  other 
hand,  goes  back  to  a  popular  incident  field  introduced  by  E.  I.  Thorsos  [28]  who — for  the 
2-D  case  derived  it  as  an  approximation  to  a  summation  of  plane  waves,  accurate  for 
sufficiently  small  angles  of  incidence  6X  (also  employed  in  [29],  [30],  [15],  [31],  [32],  [33].) 
The  limitations  of  the  2-D  scalar  Thorsos  wave  at  low  grazing  angles  were  analyzed  and 
discussed  in  [34],  [35].  The  bound  in  the  resolvability  criterion  discussed  by  Ngo  and  Rino 
[36]  also  becomes  significant  at  low  grazing  angles.  The  recommendation  for  the  3-D  vector 
case  is  to  start  over  and  simply  use  the  spectrum  given  by  (13)  which  has  the  additional 
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benefit  of  being  given  in  closed  form.  Taking  advantage  of  the  functional  dependence  of 
the  Gaussian  spectrum,  an  option  in  the  2-D  case  is  to  use  path  deformation  techniques 
to  speed  up  the  evaluation  of  the  exact  expression  for  the  incident  field  [37].  In  the  3-D 
case  we  can  at  least  band-limit  the  integration  to  a  disk  about  kp  —  kip  within  which 
the  spectrum  exhibits  a  significant  magnitude  (disk  radius  a  few  multiples  of  2 /g).  This 
leads  to  an  approximation  of  the  original  incident  field  which  satisfies  Maxwell’s  equations 
exactly.  The  derivation  of  approximate  non-Maxwellian  3-D  tapered  waves  which  can  be 
evaluated  without  integration  is  discussed  in  Section  VI. 

IV.  Polarization 


In  order  to  construct  a  wave  that  is  both  reliably  tapered  and  clearly  polarized,  for 
all  angles  of  incidence,  we  suggest  choosing  the  polarization  of  the  individual  plane  wave 
components  as  follows: 


^-hikp)  —  ■  h{kp") 

(14) 

£v{kp)  =  •  v{kp) 

(15) 

with  the  polarization  vector  of  the  central  plane  wave 

—  c^h'ip}  —  Eh  h(kip)  *f  v{kip) 

(16) 

Hence,  in  dyadic  notation 

e(kp )  =  ei  •  h(kp)  h{kp)  +  v(kp)  v(kp) j 

"  (17) 

and 

h(kp)  =  e;  •  [u(£p)  h(kp)  -  h(kp)  v(kp)} 

(18) 

The  dominant  polarization  state  of  the  tapered  wave  is  then  determined  by  the  choice 
of  Eh  and  Ev  in  (16)  which  describe  the  (in  general  elliptical)  polarization  of  the  central 
plane  wave. 

Note  that  with  this  choice  the  integrands  of  (1),  (2)  are  continuous  at  kp  =  0  [as  follows 
from  Section  VI,  we  have  in  fact  analyticity  throughout  the  kx-ky  plane  excluding  the  circle 
|A;p|  =  k  provided  an  analytic  spectrum  such  as  (13)  is  used;  at  |fcp|  =  k  the  integrands 
are  still  continuous]  as  opposed  to  the  tapered  wave  in  [1],  [2],  [4],  [9],  [8].  The  latter  wave 
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is  characterized  by  the  choice  eh(kp)  =  Eh  and  ev(kp)  =  Ev,  leading  to  rapidly  varying 
polarization  vectors  e  and  h  near  kp  =  0.  For  the  near  normal  incidence  case  this  will 
violate  the  basic  assumption  of  the  footprint  design  technique  described  in  Section  III. 
When  examined  numerically  it  is  found  that  the  approximation  of  a  prescribed,  e.g., 
Gaussian,  footprint  is  poor;  the  result  for  normal  incidence  shows  the  largest  intensity 
along  a  ring  in  the  x-y  plane  rather  than  at  the  center  (Fig.  2).  This  effect  is  also  evident 
from  the  following  consideration:  For  a  spectrum  that  satisfies  ip(kp)  —  ip(—kp)  it  can 
be  shown  that,  for  Ev  =  0,  we  have  Ei(p,  z)  =  —Ei(—p,  z)  with  the  consequence  Ei(p  = 
0,  z)  =  0  for  all  2  [Fig.  2(a)].  Similarly,  for  Eh  =  0  it  is  found  that  z  x  Ei(jp,z)  =  -ix 
Ei(-p,  z)  and  ix£i(p  =  0,  z)  =  0  [Fig.  2(b)].  Other  problems  are  leakage  of  the  intensity 
to  larger  radii  than  expected  (Fig.  2)  and  the  non-existence  of  a  clear  polarization  of  the 
wave.  By  using  (17),  (18)  these  problems  are  removed  (Fig.  3).  [The  101  x  101  tapered 
wave  field  values  for  the  results  in  Fig.  2,  3,  6,  and  7  were  calculated  using  a  summation 
of  128  x  128  plane  waves  with  a  2-D  DFT  sampling  of  the  kp  space.  The  spectrum  after 
[1],  [2],  [4],  [9],  [8]  was  calculated  using  a  2-D  FFT  algorithm.  The  horizontal  periodicity 
of  the  fields  in  the  space  domain  was  in  all  cases  30  A,  i.  e.,  twice  the  surface  length  shown 
in  the  figures,  in  order  to  avoid  aliasing  (Section  V).] 

The  tapered  wave  with  polarization  vectors  (17),  (18)  is  optimal  in  a  least  squared  error 
sense.  Consider  a  vector  field 

OO 

£i(r)  =  et  I  dkp^~p-^zH(kP)  (19) 

—  OO 


obtained  by  multiplying  a  scalar  tapered  wave  with  the  constant  polarization  vector  e, 
as  in  (16).  This  field  combines  the  desirable  properties  of  well  defined  polarization  and 
controllable  tapering.  (Note  that  | £i(p,  z  =  0)|/|ef[  corresponds  to  the  prescribed  footprint 
as  discussed  in  Section  III  and  illustrated  in  Fig.  1.)  However,  the  field  defined  by  (19) 
is  not  a  valid  electric  field  because  in  general  V  •  Si  ^  0.  We  can  therefore  ask  for  a 
permissible  wave  of  form  (1)  with  the  same  spectrum  ip  that  approximates  £i  as  close  as 
possible.  Defining 


OO 

S(z)  =  J  dp  |E,05,  z) z) 


(20) 


—  OO 
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we  find  from  Parseval’s  theorem  for  2-D  Fourier  transforms 


OO 

S(z)  =  4n2  1  dkp  e~lkzZ  ip(kp) |  | e(kp)  -e<| 

-OO 

(21) 

To  minimize  S  for  all  z  note  that 

|e(fcp)  —  e<  =| e(kp)-ei-  h(kp)h(kp)  +  v(kp)v(kp)]\  +  ei-k(kp)\ 

(22) 

where 

^  hihp)  —  ~j~  (*£  +  y  fcy  Z  ^2) 

(23) 

Hence, 

(17)  yields  the  optimal  and  the  minimum  Sm  is  given  by 

00 

Sm(z)  =  47 r2  J  dkp  e~ikzZ\p(kp)  f  \e{  •  k(kp)  2 

(24) 

—  OO 


It  is  emphasized  that  we  refer  to  S(z)  and  Sm(z)  as  “errors”  only  in  the  familiar  mathemat¬ 
ical  sense.  The  purpose  of  comparing  with  the  non-Maxwellian  field  £i(r)  is  to  uniquely 
identify  a  functional  dependence  of  e(kp)  which  can  be  expected  to  guarantee  tapering 
and  a  dominant  polarization  state  of  the  total  field  (both  as  prescribed) .  In  other  words, 
£i(r )  which  is  ideal  with  respect  to  tapering  and  polarization  is  projected  into  the  space 
of  waves  constructed  as  2-D  superpositions  of  plane  waves,  lending  its  desirable  properties 
to  an  exact  solution  of  Maxwell’s  equations. 

To  illustrate  the  approximation  behavior  numerically  we  computed  the  relative  root 
mean  squared  (RMS)  error 


where  £*  is  formed  using  the  spectrum  (13),  for  varying  tapering  parameter  g  and  incidence 
angle  6X  (in  Fig.  4  and  5  contour  levels  decrease  monotonically  for  fixed  9X  and  increasing 
g  and  are  separated  by  2dB  steps.)  The  results  in  Fig.  4  for  the  tapered  wave  in  [1],  [2], 
[4],  [9],  [8]  exhibit  the  previously  mentioned  problems  near  normal  and  grazing  incidence. 
It  is  noted  that  for  intermediate  angles  8t  and  larger  g  the  error  can  be  smaller  than 
1%  (— 20  dB)  and  that  the  approximation  behavior  for  horizontally  polarized  [Fig.  4(a)] 
and  vertically  polarized  [Fig.  4(b)]  plane  wave  components  is  similar.  For  the  tapered 
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wave  composed  according  to  (17)  and  (13)  and  for  horizontal  polarization  [Fig.  5(a)]  the 
error  is  small  everywhere  and  practically  independent  of  0t.  For  vertical  polarization 
[Fig.  5(b)]  the  error  grows  larger  towards  grazing  but  does  not  exceed  moderate  levels. 
The  fact  that  approximating  a  vertically  polarized  plane  wave  near  grazing  incidence  is 
harder  can  be  understood  intuitively  by  noting  that  the  energy  flow  of  the  tapered  wave  has 
to  “bend  down”  in  order  to  form  the  exponentially  space-limited  footprint,  a  requirement 
in  contradiction  with  maintaining  a  vertical  polarization  state.  However,  Fig.  5(b)  shows 
that  the  optimal  approximation  finds  a  reasonable  compromise.  [For  the  results  shown 
in  Fig.  4  and  5  the  RMS  error  was  evaluated  using  a  Gauss- Legendre  quadrature  over  a 
surface  of  size  7  g  x  7  g,  choosing  in  both  dimensions  5  times  the  number  of  sampling  points 
obtained  when  rounding  7  gj A  to  the  nearest  integer.  128  x  128  plane  waves  were  summed 
to  space  domain  fields  with  horizontal  periodicity  of  7  g.  The  tapering  parameter  g  was 
changed  in  steps  of  A/2  and  the  angle  of  incidence  #*  in  steps  of  5°.] 

Another  important  property  of  the  wave  based  on  (17)  is  found  from  (18)  by  noting 
that 

e*  •  h{kp)  =  h(kp)  •ei=ei-  \v(kp)  h(kp )  -  h(kp )  •O(fcp)]  •  e*  —  0  (26) 

and  thus,  according  to  (2),  e,  •  77* (f )  =  0.  The  total  magnetic  field  of  the  tapered  wave  is 
everywhere  perpendicular  to  the  electric  field  of  the  central  plane  wave. 

It  is  remarked  that  this  is  reminiscent  of  the  tapered  wave  (given  for  the  case  of  normal 
incidence  only)  in  [23],  designed  to  have  no  y  component  of  the  magnetic  field.  Setting 
e*  =  —y  in  (17),  (18)  or,  more  conveniently,  in  (34),  (37)  of  Section  VI  and  comparing 
to  (7),  (8)  in  [23]  it  is  found  that  the  waves  are  different.  In  particular,  the  polarization 
vectors  in  [23]  are  unbounded  as  — >  k  while  being  analytic  throughout  the  kx-ky  plane 

excluding  the  circle  \kp\  =  k. 

The  tapered  wave  given  previously  by  Tran  and  Maradudin  [38]  and  for  the  case  of 
vertical  polarization  employed  in  [3],  [5],  [6] — when  generalized  to  arbitrary  azimuthal 
angle  of  incidence  and  cast  into  our  formalism — turns  out  to  be  somewhat  related.  Their 
magnetic  polarization  vector  for  horizontal  polarization  is  collinear  to  (18)  when  Ev  =  0. 
However,  it  is  normalized  to  unit  length  and  the  magnetic  polarization  vector  for  vertical 
polarization  is  then  obtained  by  taking  the  vector  product  with  k(kp).  It  is  seen  that  this 
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construction  will  not  lead  to  an  optimal  approximation  of  (19)  and  thus  to  a  different 
wave. 

Finally,  we  point  out  that  our  tapered  wave  has  been  derived  by  optimizing  the  electric 
field  with  respect  to  an  ideal  field  £*.  The  magnetic  field  of  the  tapered  wave  then  followed 
from  the  familiar  relation  between  the  electric  and  magnetic  field  of  a  plane  wave  (Faraday’s 
law).  It  is  clear  that  in  a  similar  manner  we  could  derive  a  dual  tapered  wave  which  is 
obtained  by  choosing  the  magnetic  polarization  vectors  h  with  respect  to  a  non-Maxwellian 
field  Hi  and  applying  Ampere’s  law  to  find  the  electric  field. 

V.  Issues  in  the  Application  to  3-D  Scattering 

The  tapered  wave  introduced  in  this  paper  can  be  used  for  the  simulation  of  scattering 
from  randomly  rough  surfaces  with  a  planar  mean  surface.  In  a  more  complex  scenario, 
objects  are  embedded  in  a  layered  background  with  rough  interfaces.  If  the  objects  are  at 
least  partially  situated  in  the  half  space  where  the  sources  of  the  incident  wave  reside  it  is 
important  also  to  pay  attention  to  the  distribution  of  the  tapered  wave  for  z  >  0. 

Fig.  6  and  7  illustrate  the  cases  of  oblique  and  grazing  incidence,  respectively.  Fig.  6(a) 
shows  how  the  tapered  wave  forms  a  slightly  converging  beam,  approximating  the  pre¬ 
scribed  footprint  at  z  =  0  [Fig.  6(b)].  The  non- zero  intensity  in  the  top-right  corner  of 
Fig.  6(a)  is  due  to  the  periodic  nature  of  the  discretized  versions  of  (1),  (2)  with  respect  to 
p.  This  aliasing  effect,  which  in  the  present  case  would  have  no  effect  on  the  illumination 
of  objects-  situated  relatively  close  to  the  surface  at  z  =  0,  can  be  reduced — as  usual — by 
sampling  finer  with  respect  to  kp.  For  footprints  where  ip(kp)  is  not  given  in  closed  form 
as  in  (13)  but  is  computed  by  2-D  FFT  this  is  achieved  by  applying  zero  padding  before 
carrying  out  the  transformation.  The  remarkable  fact  about  Fig.  7  is  that  the  inclusion 
and  correct  treatment  of  evanescent  waves  enables  synthesis  of  the  prescribed  footprint 
even  for  9i  =  90°  [Fig.  7(b)].  Aliasing  for  2  >  0  in  this  case  is  more  severe  [Fig.  7(a)]. 

In  typical  applications  of  the  tapered  wave  concept,  electromagnetic  wave  scattering 
from  a  conducting  object  over  a  conducting  rough  surface  is  simulated  and  Glisson’s  over¬ 
lapping  triangular  flat  vector  basis  functions  [39],  [40],  [41]  for  the  electric  surface  current 
on  both  object  and  rough  surface  are  used  in  discretizing  the  electric  field  integral  equa¬ 
tion,  applying  a  Galerkin-type  method  of  moments.  We  compared  the  results  of  such  a 
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scattering  code  with  those  obtained  uy  the  hybrid  method  described  in  [10],  [11],  [12]. 
The  major  advantage  of  this  hybrid  method  is  that  the  decomposition  into  flat  surface 
problems  with  impressed  equivalent  sources  that  are  determined  by  lower  order  solutions 
allows  introduction  of  the  tensor  Green  function  for  layered  media.  This  removes  the  need 
to  solve  for  the  surface  currents  on  the  rough  surface  and  to  truncate  its  physical  dimen¬ 
sions.  In  the  comparison,  the  same  rough  surface  profile,  the  same  patch  model  for  the 
object,  and  the  same  tapered  incident  wave  were  used  for  solving  the  problem  with  the  two 
independent  codes;  reasonable  agreement  was  obtained.  Discrepancies,  however,  occurred 
for  near-grazing  angles  <9,  =  80 . . .  90°  where  the  pure  MoM  results  suffer  from  edge  effects 
due  to  the  truncation  of  the  rough  surface.  While  the  incident  wave  can  be  tapered  to  fall 
off  exponentially  towards  the  edges  the  scattered  fields  from  the  object  decay  only  as  1/r, 
giving  rise  to  problems  at  very  large  polar  angles  where  the  object  acts  as  a  reflector  that 
directs  energy  towards  the  edges.  This  indicates  the  increased  difficulty  of  the  low  grazing 
angle  rough  surface  scattering  problem  when  an  object  is  present.  In  the  hybrid  method, 
the  correct  behavior  of  the  monostatic  HH  return  at  6i  =  90°  where  in  the  flat  surface  case 
the  boundary  condition  on  the  perfectly  conducting  surface  forces  a  zero  is  guaranteed. 
It  should  be  remarked  that,  although  the  utilization  of  the  tapered  wave  concept  for  the 
method  in  [10],  [11],  [12],  which  when  implemented  up  to  and  including  the  first  order 
yields  accurate  results  for  slightly  rough  surfaces,  is  not  imperative,  it  is  still  useful  there 
because  only  finite  rough  surface  profiles  can  be  processed. 

VI.  Approximate  3-D  Tapered  Waves 

A  clear  advantage  of  the  original  2-D  scalar  Thorsos  wave  and  a  major  reason  for  its 
popularity  is  the  avoidance  of  numerical  integrations  in  the  evaluation  of  the  incident  field. 
The  price  paid  is  the  non-Maxwellian  nature  of  the  approximation  which,  as  reported  in 
[34],  can  lead  to  anomalies  in  the  computed  results  of  simulations  that  require  evaluation 
of  the  incident  field  not  only  on  the  rough  surface  as,  e.  g.,  in  object-surface  interaction 
problems.  Also  the  breakdown  of  the  approximation  near  grazing  incidence  causes  serious 
problems  in  some  applications.  Keeping  these  limitations  in  mind  and  analogous  to  the 
derivation  of  the  Thorsos  wave,  in  the  3-D  case  with  the  spectrum  (13)  we  can  argue,  in  a 
spirit  similar  to  Laplace’s  method  for  the  asymptotic  expansion  of  integrals  [42],  [43],  that 
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for  large  g  the  main  contribution  to  the  superposition  of  plane  waves  comes  from  around 
kp  =  kip.  Using  the  truncated  bivariate  Taylor  expansion  in  kx,  ky 


where  k{z  =  k  cos  we  can  obtain  from  (12)  by  carrying  out  the  integrations  and  sym¬ 
metrizing  the  result  with  respect  to  x  and  y  (without  any  further  approximations) 


e,  eiki'r  t- r  exp 
u{z) 


sir) 

g2v?(z) 


(28) 


where  k{  =  kip  —  z  kiz  and 


In  deriving  (27)-(30),  the  dispersion  relation  k2z  =  k2  —  k2x  —  k2y  has  been  used.  As 
expected,  (28)  coincides  with  (12)  exactly  when  z  =  0.  As  g  — »  oo  the  plane  wave  case 
is  recovered.  It  is  remarked  that  (28)  is  not  a  direct  generalization  of  the  Thorsos  wave 
to  3-D  because  of  the  different  formulation  of  the  superposition  integral  used  as  starting 
point.  However,  using  the  following  argument  we  arrive  at  a  condition  for  the  validity 
of  (28)  that  is  similar  to  the  one  given  in  the  2-D  case  [36],  [34],  [35],  in  particular  the 
dependence  on  (f  —  #t)2  near  grazing  carries  over  to  3-D:  The  radius  of  convergence  of  the 
full  Taylor  series  (27)  is  limited  to  k  -  kip  because  of  the  branch  point  of  the  square  root 
function.  Thus 

kg  (1  —  sin#i)  1  (31) 


is  required  for  (28)  to  be  an  accurate  representation  of  (12).  In  addition,  the  error  of  the 
truncated  series  (27)  is  multiplied  by  2  with  the  consequence  that  the  largest  |2|  considered 
should  be  small  relative  to  g ,  i.e.,  g  »  |2|max- 
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Approximations  for  the  3-D  vector  wave  case  are  derived  in  a  similar  fashion  by  addi¬ 
tionally  expanding  the  polarization  vectors.  Substituting 


h(kp)  h{kp)  =  ^[xxky-  (xy  +  yx)kxkv  +  yykl 


v(kp)v(kp)  =  j-^[xxkl  +  (xy  +  yx)kxky  +  yykl 


+  -|  [(£  z  +  z  x)  kx  +  (y  z  +  z  y)  ky]  +  z  z  t§ 


into  (17),  we  find 


S(fc»)  =  ^-[iHtf-kD+yylk'-kD  +  zikl 

-  (xy  +  yx)kxky  +  (xz  +  zx)kxkz  +  (y z  +  zy)kyk: 


Equation  (34)  also  follows  from  hh  +  iiv  =  I  —  kk  with  I  the  identity  tensor  and  k  as  in 
(23).  Similarly,  with  the  dyads 

v(kp)h(kp)  =  [{xx-yy)kxky-  xykl  +  yxk^ 

+  ^(zxky  -  zykx)  (35) 

h(kp)  v(kp)  =  [(s  x  -yy)kxky  +  xyky  -yxkl] 

+  ^(xzky-yzkx)  .  .  (36) 

we  obtain  from  (18) 


h(kp)  ‘  =  ■  [{xy  -  yx)kz  +  (x  z  -  zx)  ky  -  (yz  -  zy)  kx\  (37) 

It  is  observed  from  (34),  (37)  that  the  tensors  in  (17),  (18)  are  symmetric  and  anti¬ 
symmetric,  respectively.  More  importantly  here,  it  is  evident  from  (34),  (37)  that  both 
e(fcp)  and  h(kp),  viewed  as  functions  of  the  two  real  variables  kx  and  ky,  are  analytic 
throughout  the  kx-ky  plane  excluding  the  1-D  set  of  points  forming  the  circle  \kp\  =  k.  (It 
is  emphasized  that  the  region  of  analyticity  includes  kp  =  0,  c.  f.,  Section  IV.)  Thus,  the 


Taylor  series 


e(kip  +  kp)  —  ^2  amn  k”1  ky 

771,71=0 
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1  [  dm+n  x 

amn~  m\n\  dk?dk"e{  p)  -k  r 

L  x  y  J  kp=ki. 


converges  in  the  disk  |/cp|  <  k  -  kip,  and  similar  expressions  hold  for  h(kp).  Applying 
the  approximation  (27)  to  (1)  with  the  spectrum  (13),  inserting  (38),  and  using  twice  the 
Fourier  integral  identity  (c.  f.,  3.958.2  in  [44]) 


dxxneiqx~ax 


TT  dn  _sl 

i  n  —  e  4a 


a  dqn 


fW  /  i  \ 
V  a  \2aJ 


n  _a£  ^  (~fl)m 
6  4°  i'0(n-2m)!m! 


where  [n/2]  is  the  integral  part  of  n/2,  leads,  again  without  further  approximation,  to  an 
expression  of  the  form 


Ei(r)  «  e,fc<  r  — exp  — *  r  amnpmn(r) 

u{z)  [  g2v?(z) J  m%0 


where  s(r)  and  u(z)  are  as  in  (29),  (30),  the  amn  are  obtained  from  (34),  and  pmn(f)  is 
polynomial  in  x,  y,  z  up  to  ^-dependent  correction  factors  that,  similar  to  (30),  are  unity 
for  2  =  0  or  g  -*  oo.  It  is  noted  that,  different  from  the  scalar  case  discussed  above,  (41) 
for  2  =  0  is  only  an  approximation  of  the  superposition  of  plane  waves  that  we  started 
out  with  because  of  the  finite  radius  of  convergence  of  (38)  and  the  fact  that  the  spectrum 
-  (13)  is  not  band-limited.  The  conditions  on  g  for  (41)  to  be  a  reasonable  approximation 
are  as  stated  above  for  the  scalar  case.  We  have  Goo  =  e*  and  poo  —  1-  Thus  the  lowest- 
order  electric  field  approximation  following  from  (41)  is  the  same  as  what  is  obtained  by 
approximating  kz  in  (19)  with  the  help  of  (27)  or  multiplying  (28)  by  e£/e£.  The  pmn  of 
higher  order  vanish  as  g  —*  oo.  The  algebraic  details  for  the  pmn  of  any  order  (integrating 
over  ky  first)  are  as  follows: 


pmn(r)  = 


n—2k—j 


,[  2i  ]"  ^  hrtWf 

t„  4<=fc! 

x  ny  _ + _ {>/  +  — 

jfe  J-(n  -2k-j)l  V  kit 


x  Y  [  \(z)\ _ \x+hlz  +  w(y  z) 

to  4<(rn  +  j-2t)W.  [X  +  kiz  +  Hv' Z) . 


k ix  k{y 


9 2  u\{z) 


m+j— 2t 
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where 


VII.  Conclusion 

We  considered  the  problem  of  constructing  a  3-D  tapered  wave  as  a  superposition  of 
plane  waves,  taking  into  account  both  propagating  and  evanescent  waves.  The  use  of  the 
simple  Gaussian  plane  wave  spectrum  was  recommended  in  order  to  avoid  problems  near 
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the  grazing  incidence.  The  introduced  special  choice  for  the  polarization  vectors  removed 
the  problems  of  losing  a  dominant  polarization  state  and  degradation  of  tapering  near 
the  normal  incidence.  Mathematically  speaking,  the  proposed  polarization  vectors  are 
analytic  at  the  origin  of  the  2-D  wavenumber  space.  Moreover,  the  choice  of  polarization 
vectors  was  shown  to  lead  to  an  exact  solution  of  the  Maxwell  equations  which  is  an 
optimal  approximation  of  an  ideal  but  non-Maxwellian  tapered  field  that  is  constructed 
by  multiplying  a  scalar  tapered  wave  with  a  constant  polarization  vector.  The  result  is 
a  reliably  tapered  wave  with  a  dominant  polarization  state  that  can  be  used  uniformly 
for  all  angles  of  incidence.  We  discussed  the  application  of  the  proposed  tapered  wave 
in  simulating  3-D  electromagnetic  scattering  from  a  conducting  object  over  a  conducting 
rough  surface.  Newly  encountered  problems  near  the  grazing  incidence  were  attributed 
to  secondary  edge  effects  which  are  unrelated  to  the  tapered  incident  wave  but  indicate 
the  difficulty  of  the  low  grazing  angle  rough  surface  scattering  problem  when  objects  are 
present.  It  was  pointed  out  that  methods  which  avoid  such  edge  effects  could  also  benefit 
from  the  utilization  of  the  tapered  wave.  In  some  situations  it  might  be  desirable  to  have 
an  approximate  3-D  tapered  wave  at  one’s  disposal  which  does  not  require  a  2-D  numerical 
integration  (summation  of  plane  waves),  trading  accuracy  in  satisfying  Maxwell’s  equations 
for  computational  speed.  We  presented  the  derivation  of  approximations  for  both  the  3- 
D  scalar  and  vector  case.  The  expansion  of  the  polarization  vectors  is  based  on  their 
analyticity.  The  local  character  of  the  technique  employed  forces  the  breakdown  of  the 
approximations  at  grazing  incidence. 
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Fig-  1. 


Example  of  a  prescribed  Gaussian-shaped  footprint  to 


be  approximated  by  the  vector  tapered 


wave  ( g  =  2  A). 
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(b)  Vertically  polarized  plane  wave  components. 

Fig.  2.  Resulting  footprints  at  normal  incidence  for  the  tapered  wave  after  [1],  [2],  [4],  [9],  [8].  The 
approximation  of  the  prescribed  footprint  (Fig.  1)  is  not  satisfactory. 
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(b)  Vertical  polarization. 

Fig.  3.  Resulting  footprints  at  normal  incidence  for  the  tapered  wave  introduced  in  the  present  paper 
and  approximating  the  prescribed  footprint  of  Fig.  1. 
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(a)  Horizontally  polarized  components.  (b)  Vertically  polarized  components. 

Relative  RMS  error  [dB]  at  z  =  0  for  the  tapered  wave  after  [1],  [2],  [4],  [9],  [8]  as  compared  to  a 
-Maxwellian  field  with  prescribed  tapering  and  polarization. 


(a)  Horizontal  polarization. 


(b)  Vertical  polarization. 


Fig.  5.  Relative  RMS  error  [dB]  at  z  =  0  for  the  tapered  wave  introduced  in  this  paper. 
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(a)  Side  view  (x  =  0). 


(b)  Footprint. 


Fig.  6.  Beam  formation  of  the  tapered  wave  at  oblique  incidence  (8t  =  40°,  fa  =  90°,  g 
polarization). 


2  A,  horizontal 
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(a)  Side  view  (x  =  0). 


Fig. 


(b)  Footprint. 


7.  Beam  formation  of  the  tapered  wave  at  grazing  incidence  (0*  =  =  90°,  g  =  2  A,  horizontal 

polarization). 
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